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ABSTRACT 

We present new results on the kinematics and spatial distribution of metal-enriched gas within ~ 125 kpc of 
star-forming ("Lyman Break") galaxies at redshifts 2 £ z ^ 3. In particular, we focus on constraints provided 
by the rest-frame far-UV spectra of faint galaxies- and demonstrate how galaxy spectra can be used to obtain 
key spatial and spectral information more efficiently than possible with QSO sightlines. Using a sample of 89 
galaxies with (z) = 2.3 ±0.3 and with both rest-frame far-UV and Ha spectra, we re-calibrate the measurement 
of accurate galaxy systemic redshifts using only survey-quality rest-UV spectra. We use the velocity-calibrated 
sample to investigate the kinematics of the galaxy-scale outflows via the strong interstellar (IS) absorption lines 
and Lyman a emission (when present), as well as their dependence on other physical properties of the galaxies. 
We construct a sample of 512 close (1 - 15" ) angular pairs of z ~ 2-3 galaxies with redshift differences 
indicating a lack of physical association. Sightlines to the background galaxies provide new information on 
the spatial distribution of circumgalactic gas surrounding the foreground galaxies. The close pairs sample 
galactocentric impact parameters 3-125 kpc (physical) at (z) = 2.2, providing for the first time a robust map of 
cool gas as a function of galactocentric distance for a well-characterized population of galaxies. We propose 
a simple model of circumgalactic gas that simultaneously matches the kinematics, depth, and profile shape 
of IS absorption and Lya emission lines, as well as the observed variation of absorption line strength (H I 
and several metallic species) versus galactocentric impact parameter. Within the model, cool gas is distributed 
symmetrically around every galaxy, accelerating radially outward with v out (r) increasing with r (i.e., the highest 
velocities are located at the largest galactocentric distances r). The inferred radial dependence of the covering 
fraction of cool gas (which modulates the absorption line strength) is f c (r) oc r" 7 with 0.2 7 0.6 depending 
on transition. We discuss the results of the observations in the context of "cold accretion", in which cool gas 
is accreting via filamentary streams directly onto the central regions of galaxies. At present, we find little 
observational evidence for cool infalling material, while evidence supporting the large-scale effects of super- 
wind outflows is strong. This "pilot" study using faint galaxy spectra demonstrates the potential of using 
galaxies to trace baryons within galaxies, in the circumgalactic medium, and ultimately throughout the IGM. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Observational studies of the galaxy formation process are 
reaching a critical juncture, where the accumulation rate of 
new data may be overtaking our ability to reach new under- 
standing. There are now many large surveys designed to study 
galaxy formation and evolution over ever-increasing volumes 
and redshift ranges, backed by the unquestionable power of 
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multi-wavelength observations from X-rays to the far-IR/sub- 
mm and radio. Interpretation of the survey data is aided by 
a multitude of theoretical results based on simulations and/or 
semi-analytic calculations. However, many of the most funda- 
mental remaining questions in galaxy formation involve com- 
plex baryonic processes that are difficult to model, are not 
well-constrained by current observations, or involve physics 
that are not yet well understood. 

While there is general agreement about the development of 
structure in the dark matter component on scales larger than 
that of galaxies, the astrophysics of the baryonic response 
to the dark matter structure, the subsequent feedback of en- 
ergy from star formation, supernova explosions, and accre- 
tion energy from the growth of supermassive black holes, 
and the flow of gas into and out of galaxies, remain largely 
unconstrained- and thus subject to substantial debate. "Feed- 
back" has become a buzzword, universally acknowledged as 
something important to understand, but there is little agree- 
ment about what it really means and how it affects the "big 
picture". Nevertheless, some kind of feedback is invoked to 
explain many otherwise-inexplicable observations: the cessa- 
tion of star formation in massive galaxies at high redshift; the 
small number of low-mass galaxies relative to the dark mat- 
ter halo mass function predicted by the otherwise successful 
ACDM cosmology; the correlation between galaxy spheroid 
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mass and the mass of central super-massive black holes; the 
general absence of cooling flows in clusters of galaxies; the 
metal enrichment of intracluster gas, just to name a few. The 
most obvious sources of the energy and/or momentum re- 
quired to explain these phenomena are massive stars, super- 
novae, and AGN activity, all of which must have exerted most 
of their influence in the distant past, when these processes 
were at their peak intensity. 

One route to understanding the relevant baryonic processes 
is via simultaneous study of galaxies and the IGM, in the 
same cosmic volumes during the epoch when they were ar- 
guably exerting the greatest influence on one another- near 
the peak of both universal star formation and supermassive 
black hole growth in the redshift range 3 ^ z ^ 1.5. Com- 
bining two powerful lines of investigation provides comple- 
mentary information on the state of baryons, both those col- 
lapsed into galaxies and those residing outside of galaxies. 
The IGM and "circumgalactic medium" (CGM; by which we 
mean the gas-phase structures found within 300 kpc [phys- 
ical] of galaxies), together present a laboratory in which the 
effects of galaxy formation and AGN accretion (e.g., radiative 
and hydrodynamical "feedback" and its recent history) can be 
measured on scales that are not accessible using direct obser- 
vations of galaxies. Similarly, the galaxy distribution relative 
to the lines of sight to background objects tells us more about 
how the physical information garnered from the absorption 
line studies should be interpreted. These ambitious science 
goals - to observe both diffuse gas and galaxies/ AGN in the 
same survey volume at high fidelity and down to small scales 
- require a different approach as compared to most spectro- 
scopic surveys of the distant universe. Whereas the move- 
ment of most galaxy surveys has been toward larger and larger 
scales, the most vexing remaining uncertainties are related to 
phenomena occurring on ~ 1 - 10 Mpc scales, where the fig- 
ure of merit is information density and not total survey vol- 
ume. 

It has been known for some time that galactic-scale out- 
flows with velocities of several hundred km s" 1 are ubiq- 
uitous in star-forming galaxies at all redshifts for which 
interstellar (IS) ab s orption fea t ures a re accessible (e.g. , 
I Steidel et all 11996b lFranx etal.1 I l997t IPettini et all 120001. 
20011: IShaolev et alJl2003b iMartinl 120051: iRupke etal] £0051: 
iTremonti et alJ l2007: Wei ner et alJ l2009.) Evidence is found 
in the offsets between the redshifts of the nebular emis- 
sion lines, interstellar absorption lines, and Lya emis- 
sion (IPettini et alj 12001). in the relative velocities of stel- 
lar, interstellar and nebular lines in composite UV spectra 
( Shapl ey et alJ l2003). and in the correlation of C IV systems 
seen in absorption i n QSO spectra with the positions of the 
galaxies themselves ( Adelberger et al. 2003; Adelberger et al. 
2005a). The outflowing gas has been obse rved in de- 
tail for only a handful of lensed galaxies (IPettini et alj 
I2002t lOuider et al1l2009allbl) : in the best-observed of these, 
MS1512-cB58, the interstellar absorption is distributed over 
a very wide range of velocities, -800 ^ v ^, + 200, with a 
centroid velocity of -255 km s" 1 ( IPettini et alj|2002l) . Such 
outflows are also a general feature of starburst galaxies in 
the local universe, where rest-UV spectroscopic studies re- 
veal similar ly complex velocity s t ructure and gas in multiple 
phases (e.g. iHeckman et allll990l:lLehnert & HeckmarH l 1996 1: 

ISchwartz et alj [2006; 
However, the outflow 
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phenomenon is so much more widespread at high redshifts 



that it is influencing essentially every galaxy, and potentially 
every galaxy's local environment. It must have an influence 
on both the chemical evolution of the galaxies and of the IGM, 
it may well regulate the maximum star formation rate attain- 
able by a galaxy, and without doubt it is an essential ingre- 
dient to basic understanding of the circulation of baryons as 
galaxies are forming. 

In spite of the substantial observational evidence for 
galaxy-scale outflows, most of the recent theoretical work has 
focused instead on the infall of cool gas ("Cold Accretion" 
or "Cold Flows") via filaments, directly onto the central re- 
gions of forming galaxies. According to much of the recent 
work, this mode of gas accretion is what feeds (and regu- 
lates) high star formation rates in high redshift galaxies un- 
til they attain a particular mass threshold (M tot > 10 12 M ) at 
which point virial shocks develop and accretion of cold mate- 
rial is suppressed (iKeres et alJl2005ll2009t lDekel et alJl2009l: 
iBrooks et al] 12009c ICeverino et al. 2009j). This transition is 
believed to be at least partially responsible for producing mas- 
sive "red and dead" galaxies as early as z ~ 2.5. Similarly, 
the phenomenon of spati ally extended Lya emission, includ- 
ing giant Lya "blobs" ( Steidel et al. 2000) has been ascribed 
to cooling radiation from the denser portions of the cold 
streams as they a re accreting directly onto the galaxy's cen- 
tral regions (e.g., iHaiman eta l. 2000; F urlanetto et al.l l2005h 
iDiikstra & Loebl2009UGoerdt et al.ll2009l) . Although the cold 
accretion picture may be attractive, the predictions of the 
observational consequences for CGM gas are rather model- 
dependent for both absorption lines and Lya emission. The 
simulations which predict the cold accretion generally do not 
account for IS gas that may have been carried to large galac- 
tocentric radii by outflows, nor for Lya photons produced by 
photoionization in a galaxy's H II regions which then scatter 
their way in space through CGM gas before escaping. It may 
well be that the observational signatures of infall by cold ac- 
cretion and outflows via supernova-driven winds are very sub- 
tle, and perhaps indistinguishable. Possibly the most telling 
differences would be kinematic- absorption line signatures 
of infalling material in galaxy spectra would be expected to 
be primarily redshifted, while outflowing material would be 
blue-shifted and could reach much higher velocities with re- 
spect to the galaxy systemic redshift. Thus, accurate determi- 
nation of the galaxy systemic redshift is an essential part of 
understanding the nature of the CGM. 

The spatial distribution of the blue-shifted high-velocity 
material seen in absorption against the galaxy continuum is 
not yet established. We know that virtually every z > 2 galaxy 
bright enough to be observed spectroscopically is driving out 
material at velocities of at least several hundred km s" 1 , but 
we do not know how far this material travels, or even where 
it is with respect to the galaxy as we observe it. The mass 
flux associated with such flows has bee n measured in onl y 
one case at high redshift, MS1512-cB58 ( IPettini et al.ll2000l) . 
and even in this case the result depends very sensitively on the 
assumed physical location of the absorbing material responsi- 
ble for the bulk of the observed absorption. The best hope for 
constraining the location of outflowing gas is by observing ob- 
jects lying in the background, but at small angular separation, 
relative to the galaxy of interest. In this case the challenge 
is to find background objects bright enough in the rest-frame 
far-UV but close enough to the foreground galaxy to provide 
interesting constraints. A great deal of effort, over a large 
range of redshifts, has been invested using QSOs as back- 
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ground sources, where absorption by metallic ions or HI in the 
spectrum of the QSO is compared with g alaxies with known 
redsh ifts and projected separations ("e.g.. iBergeron & Boisse 
[l99TtlSteidel et alJl994tlChen et dJ2001tlSteidel et al.l2002t 
Lanzetta et alj 119951; iDanforth & Shulll 20081; iBowen et alj 
ll995tlBouche et al.ll2007l;lKacprzak et alj|2009l) : most of this 
work has focused on redshifts z < 1 because of the increas- 
ing difficulty obtaining spectra of the foreground galaxies at 
higher redshifts. Even if galaxies are identified and have red- 
shifts that correspond closely with observed absorption, the 
association of particular absorption systems with identified 
galaxies is almost always ambiguous, since the dynamic range 
for identifying faint galaxies is limited, and often the bright 
background QSOs make it challenging to observe galaxies 
within a few arcsec of the QSO sightline. Finally, there is the 
controversial issue of whether metals seen near to, but outside 
of, galaxies are a direct result of recent star formation or AGN 
activity, or are simply tracing out regions of the universe in 
which some galaxies, perhaps in the distant past, polluted the 
gas with metals (e.g.. iMadauet alJl200U IScannapieco et al. 
2002 1 ; iMori et alj 12002 iFerraral 120031 ; cf. lAdelberger et al. 
20031; lAdelberger et alj|2005al) . 

In this paper, our goal is to try to understand the spatial 
distribution of cool gas seen in absorption against the stel- 
lar continuum of every galaxy observed at high redshift. The 
objective is to empirically track the kinematics and structure 
in the CGM from the central parts of galaxies all the way to 
large galactocentric radii. In this work, we use only galaxy 
spectra, primarily in the rest-frame far-UV, but we calibrate 
the velocity zero-point using a set of nearly 100 Ha mea- 
surements in the observed-frame near-IR for galaxies in the 
redshift range 1.9 £z 2.6. The ne ar-IR measureme nts are 
drawn primarily from the sample of Er b et al.1 (|2006c), after 
which we use a much more extensive set of rest-UV spectra 
from a nearly completed UV-selected galaxy survey targeting 
the same range of redshifts. 

The paper is organized as follows: in $2] we examine 
the statistics of the kinematics of the outflows using a sub- 
sample of galaxies for which both near-IR nebular Ha spec- 
troscopy and reasonably high-quality optical (rest-UV) spec- 
tra are available. We also present new empirical formulae for 
estimating galaxy systemic redshifts for the typical case in 
which only low S/N rest-UV spectra are available. In $3] we 
seek correlations between the interstellar absorption line kine- 
matics, in particular the bulk velocities measured from strong 
low-ionization transitions, and other measured galaxy proper- 
ties, ^describes further inferences on the structure and kine- 
matics of outflowing material from high S/N composite far- 
UV spectra, while ^examines the observed behavior of Lya 
emission and its relationship to the IS absorption features, and 
attempts to understand the observations with simple models. 
We introduce in §|6]the use of close angular pairs of galaxies at 
different redshifts for mapping the spatial distribution of the 
circumgalactic gas around the foreground galaxies, while f7] 
develops a simple geometrical and kinematic model for out- 
flows consistent with both the line profiles in galaxy spectra 
and the larger-scale distribution of gas in the CGM. $8] dis- 
cusses the observational results and their interpretation in the 
context of the models, and ^summarizes the conclusions and 
discusses the prospects for improvement in the future. 

We assume a A-CDM cosmology with Cl m = 0.3, ft\ = 0.7, 
and h = 0.7 throughout, unless specified otherwise. 



2. BULK OUTFLOW VELOCITIES ASSOCIATED WITH Z ~ 2 - 3 
GALAXIES 

2.1. Spectra in the Ha Galaxy Sample 

The rest-frame far-UV spectra of star forming galaxies in- 
clude numerous absorption features which (in principle) pro- 
vide detailed information on the young OB stars responsi- 
ble for the bright continuum, as well as the IS and circum- 
galactic a tomic and ionized gas a ssociated with the galaxy 
(see, e.g., lPettini et al.ll2000l [20021) . The strongest stellar fea- 
tures are due to stellar winds from massive stars, producing 
broad P-Cygni profiles in higher ionization lines such as N V 
AA1238,1242, C IV AA1548,1550, and Si IV AA1393,1402. 
The absorption is broad (Av £ 2000 km s" 1 ) and shallow, 
with depth dependent on the metallicity of the O-stars. Pho- 
tospheric lines from the same OB stars are also present, but 
the lines are generally much weaker than the wind features, 
and so can be difficult to discern in the spectra of individ- 
ual high redshift galaxies. The IS lines are superposed on 
the integrated stellar spectrum, where absorption features of 
abundant species (e.g., H I, O I, C II, C IV, Si II, Si III, Si IV, 
Fe II, Al II) can be extremely strong — strong enough to be 
useful for redshift identification in low S/N spectra. Unfor- 
tunately, the most accessible absorption features- the IS and 
wind features - are not useful for accurate measurements of 
galaxy redshifts because of the non-gravitational origin of 
their kinematics. Similarly, Lya emission, which is observed 
in ~ 50% of a con t inuum-selected sample at z ~ 2 - 3 (e .g. 
iShaplev etal.ll200l ISteidel et alJl2004t iKornei et alJl2009t) is 
resonantly scattered, altering the emergent kinematics in a 
manner that depends on both the optical depth and velocity 
distribution of the gas doing the scattering. 

Fortunately, nebular emission lines lines originating in a 
galaxy's H II regions are also relatively accessible observa- 
tionally, with the strongest lines (e.g., Ha , H/3 , [OIII]) 
found in the observed-frame near-IR at the redshifts of in- 
terest. The Ha line, which is not resonantly scattered and 
whose strength is strongly dependent on both the ionizing UV 
radiation field and the density, is likely to provide a reason- 
able estimate of the systemic redshift of the stars, as well as a 
measure of gravita t ionally-induced motion within the galaxy 
dPetfini et al]|2001t lErb et al. 2006c]). The only disadvantage 
of the nebular spectra is that their numbers are currently small 
in comparison with the available rest-UV spectra. 

The quality of rest-UV survey spectra for typical high red- 
shift galaxies is generally not sufficient for detailed spectral 
analysis, and so much of what we know about general trends 
between spectral properties and other physical parameters 
of galaxies is based on composites ( e.g., ISteidel et al.l 1200 it 
Shaplev et al. 2003; Erb et al. 2006a). In order to form these 
composites for the far-UV spectra, one must assume a rela- 
tionship between the velocities of measured spectral features 
and the object's systemic redshift z sys ; as discussed above, this 
task is made difficult by the fact that the strong lines of reso- 
nance transitions in the far-UV are seldom closer than a few 
hundred km s" 1 to the objects' true redshifts. 

With the current, enlarged sample of z ~ 2 galaxies, includ- 
ing generally higher-quality UV spectra than available in the 
past, we revisit the measurement of galaxy systemic redshifts 
using far-UV spectral features , building o n previous results 
(Pettin i et all 1200 It lAdelberger et al.1 120031: lAdelberger et al.l 
2005a). In this section, we examine in detail the relation- 
ship between observed far-UV IS lines, Lya emission, and 
the redshifts defined by the Ha line in the rest-frame optical, 
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including careful attention to possible systematic errors. 

The redshifts (both Ha and UV-based) for all b ut 8 of the 
galaxies have been presented in Er b et al.l d2006cl) . where the 
observations are described in detail. The new H a spectra were 
obtain ed in the same manner using NIRSPEC dMcLean et al.l 
119981) in 2005 June. Using a 0'/76 entrance slit, the typical 
spectral resolution achieved with NIRSPEC is R ~ 1400, or 
FWHM~ 215 km s" 1 , so that the typical Ha line widths of 
<7 V ~ 100 km s" 1 (FWHM~ 240 km s" 1 ) are marginally re- 
solved. The rest-UV spectra were obtained with the LRIS-B 
spectrograph on the Keck 1 telescope, primarily using a 400 
line/mm grism blazed at 3400 A with a dichroic beam splitter 
sending all wavelengths A < 6 800 A into the blue channel, as 
described in ISteidel et al.l (120041) . As discussed in detail be- 
low, the typical effective spectral resolution achieved in this 
mode is ~ 370 km s" 1 (FWHM), or a res ~ 160 km s" 1 . 

For objects whose UV spectra show detectable interstellar 
absorption lines, we have included only those measurements 
deemed of sufficient quality to yield reasonably precise mea- 
sur ements. We have also excluded the handful of objects from 
the lErb et al.l l2006d sample with z < 1.9, in order to avoid 
systematic biases caused by differential atmospheric refrac- 
tion (discussed below), and to focus on the redshift range for 
which we have the highest quality ancillary data. 

The resulting Ha sample consists of a total of 89 galax- 
ies with (z) =2.27 ±0.16, of which 48 (54%) have UV spec- 
tra in which only an interstellar absorption line redshift was 
measured, 38 (43%) have measurements of both Lya emis- 
sion and interstellar absorption, and 3 (3.5%) have only Lya 
emission redshifts measured. These fractions are compara- 
ble to those obtained from the full spectroscopic sample of 
~ 1 600 galaxies in the same range of redshifts, 1.9 < z < 2.6 
(see ISteidel et al.ll2004 . Most (87% ) of the galaxies in the 
parent sample satisfy the "BX" photometri c selection criteria 
( Adelber g"eretani2004l:[Steidel et al.ll2004l) . with the remain- 
der sat isfying the "MD" color criteria defined in ISteidel et alj 
(2003). These galaxies selected on the basis of their UV color 
have been shown to include all but the dustiest star-forming 
galaxies in the redshift ra nge of interest ( Adelberger & Steidell 
l2000llReddv etani2008h . 

2.2. Redshift Uncertainties 

As discussed in detail by Shapley et al. (2003), a large num- 
ber of interstellar absorption lines is commonly observed in 
far-UV spectra of star-forming galaxies. In general, all of 
the detected lines are used to verify that the assigned redshift 
is correct, but the positions of only 3 lines: C II A 1334.53, 
Si IV A1393.76, and Si II A1526.72, have been used to mea- 
sure zis- Of the strong resonance lines in the rest wavelength 
range 1000-1600 A (the range in common for almost all of 
the spectra discussed in this paper), these lines are the least 
likely to be blended with other strong lines and most likely 
to yield a consistent measure of zis- Clearly, the accuracy of 
each measurement of zis depends on the strength and width 
of the lines, and on the quality of the spectra. The S/N of the 
spectra used for the current study varies considerably. Empir- 
ically, from internal agreement of different absorption lines 
in the same spectrum, and from repeated measurements of 
the same galaxy on different slitmasks, we estimate that the 
typical measurement uncertainties are ~ 100 km s" 1 for zis 
redshifts and ~ 50 km s" 1 for zi_ ya - 

Unfortunately, systematic uncertainties can be more diffi- 
cult to quantify. For example, the measured wavelengths of 



features depend on the illumination of the spectrograph slit 
by the object, whereas the wavelength solutions are deter- 
mined from calibration lamps and night sky spectra which 
illuminate the l'/2 slits uniformly. This source of error can 
also be wavelength-dependent due to the effects of differ- 
ential atmospheric refraction. The features used for redshift 
measurement generally fall in the observed wavelength range 
~ 3750 - 5200 A for the Ha sample considered here, with 
(zna) = 2.27 ±0.16. Over this wavelength range, the differ- 
ential refraction 13 would be ~ 0.7 arcsec at an airmass of 
secz= 1.3. Slitmasks were almost always observed within 
~ 10-15 degrees of the parallactic angle when the airmass 
was significantly different from 1 .0, but even so, the ampli- 
tude of differential refraction perpendicular to the slit could 
be as large as ~ 0.2 arcsec, which could map into a velocity 
shift of up to ~ ±75 km s" 1 depending on the seeing and the 
spatial profile of the galaxy. 

The measurements of ZHa are subject to a different set 
of systematic uncertainties. Differential refraction in the 
near-IR is negligible for our purposes, but the observations 
were obtained using a 0'.'76 slit after applying a blind off- 
set from a nearby star to the position of the galaxy mea- 
sure d from either an optic al or near-IR continuum image 
(see lErb et all I2003L l20"06ch . While the observations were 
taken in a way that should minimize any systematic offsets 
in the resulting zn a due to pointing or astrometric inaccura- 
cies, they of course measure the velocity of only the flux that 
entered the slit. Empirically, we found that repeated obser- 
vations of the same galaxy (using different slit position an- 
gles) suggest an accuracy of ±60 km s" 1 for zho (nns), with 
the largest excursions from consistency applying to objects 
known to be spatially extended. Another direct comparison 
of the redshifts for 14 objects from the NIRPSEC sample of 
lErb et al] d2006d) is provided by observations of the same ob- 
j ects with the integral-field spe ctrometer SINFONI at the VLT 
dForster Schreiber et alj 12006). The SINFONI redshifts are 
based on the velocity centroid of all detected Ha flux from the 
object, and are not affected significantly by pointing errors or 
by slit losses. This comparison shows a level of agreement 
similar to our estimate from multiple observations with NIR- 
SPEC: (zn-Zs) = -34 ±59 km s _1 (rms), where zn is the NIR- 
SPEC redshift and z s is the SINFONI redshift. The average 
offset is marginally significant (~ 2a), but the scatter is con- 
sistent with our NIRSPEC experience. Eleven of the galaxies 
in the present Ha sample were also observed using OSIRIS ( 
with Laser Guide Star Adapti ve Optics (LGS AO) on the Keck 
2 telescope dLaw et alj|2009l) : the average redshift difference 
Zn~Zo = 8 ±41 km s" 1 , indicating no evidence for a system- 
atic difference. 

In summary, adopting an uncertainty in zho of 60 km s" 1 , 
the typical uncertainties in the measured values of Avis = 
c (z IS - z Hq ) /( 1 + z Hq ) and A v Lyce = c(z Lya - z Ha )/ ( 1 + z Ha ) f or 
a given galaxy are ~ 130 km s" 1 and ~ 90 km s" 1 , respec- 
tively. 

2.3. Sample Statistics 

Fig.Q]shows a histogram of the IS and Lya emission veloc- 
ities with respect to the systemic redshift defined by ZHa- The 

13 All of the spectroscopic observations of the Ha sample were obtained 
before the commissioning of the Cassegrain ADC on Keck 1. We are ex- 
cluding for the moment galaxies at redshifts where key features fall at wave- 
lengths shorter than 3750 A due to the rapidly increasing amplitude of differ- 
ential refraction. 
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FIG. 1. — Histogram of the measured (centroid) velocities of interstellar 
absorption lines (blue) and Lyman a emission (red) with respect to the galaxy 
nebular redshift as defined by the centroid of the Ha emission line, for a 
sample of 89 galaxies with (z) = 2.27 ± 0. 16. The sample includes only those 
galaxies having both nebular line redshifts and rest-UV spectra of adequate 
quality to measure absorption line centroids. In this sample, 86 of the 89 
galaxies have measured values of zis, 3 have only Zhya, and 39 have both. 
The mean values of the velocity offsets are indicated. 

distributions have (Avis) = — 1 64 ±16 km s" 1 , and (AvL yQ ) = 
445 ± 27 km s" 1 , where the quoted uncertainties are errors 
in the mean. Fig. [2] illustrates the same sample as in Fig.Q] 
where different symbols are used depending on the UV spec- 
tral morphology of the galaxies. 

In light of the current Ha sample of z ~ 2.3 galaxies, it 
is worth re-examining the "rules" that one would use to es- 
timate the true systemic redshift of the galaxies given only 
information contained in their rest-UV spectra, and assuming 
that zho defines the res t frame. Using a l i near reg ression form 
similar to that used by Adelberger et al. (2005a), for galaxies 
with both Zhya and zis measurements, 

z Ha =Zis + 0.00289- 0.0026(2.7 -z IS ) a z = 0.00127 (1) 

ZHa=ZLya-0.0054+0.0001(2.7-z Lya ) <7- = 0.00193, (2) 

corresponding to velocity offsets of Avrs = -170 ± 115 
km s" 1 and Avl yq = +485 ± 175 km s" 1 , respectively, at the 
mean redshift of zn a = 2.27. For objects with a measurement 
of zis only, 

ZHa=zis+0.00303-0.0031(2.7-zis); a 2 = 0.00145, (3) 

or Avis = — 165 ± 140 km s" 1 (error is the standard deviation) 
at the mean redshift of the sample. 

Using all 86 Ha galaxies with measured zis, the best fit 
single relation of the form in eqs. 1-3 is 

Zito=zis + 0.00299-0.00291(2.7-z IS ) <r z = 0.00138 (4) 

or Avis = -166 ± 125 km s" 1 at (zhc) = 2.21. We find that 
including z LyQ in the above regression formulae increases the 
rms redshift uncertainty over that obtained using only zis, in 
contrast to similar estim ates at somewhat higher redshift by 
lAdelbe rger et al. (2005a). One possible explanation for the 
difference could be the generally weaker Lya lines in the 



z ~ 2.3 sample compared to that at z — 3 dReddv et alj '2008). 
We return in ^5] to a discussion of the kinematics of the Lya 
emission line. In any case, using only the absorption redshift, 
with a constant offset of ~ +165 km s -1 , would provide an 
estimate of zna accurate to ~ 125 km s" 1 (rms). 

There are too few objects (3 of 89) in the Ha sample hav- 
ing only Zhya to define a significant relationship for such ob- 
jects (which are also quite rare in the full z ^ 2 spectroscopic 
sample), although these 3 objects have (Avl Y q) = 400 ± 183, 
consistent with equation [2] above. For this reason, we use 
equation [2] for subsequent estimates of zua when only Lya 
emission is available. 
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Velocity Relative to z(Ha) 

FIG. 2. — Plot showing the interstellar absorption line centroid velocities 
(blue) and centroid Lyman a emission velocities relative to the redshift de- 
fined by Ha for the same sample as in Fig. [T] Galaxies for which both in- 
terstellar absorption redshifts and Lyman a emission redshifts are available 
are indicated with blue (absorption) and red (emission) solid dots; open trian- 
gles show systems for which one or the other measurements is not available. 
The circle d objects are one s which exhibited measurable velocity shear in the 
sample of Erb et al. 2006c (see text for discussion). 

Figs. [T]and [2] show that a significant fraction of the galax- 
ies have IS absorption line centroid velocity shift Avis con- 
sistent with zero. Given the uncertainties in Avis, this is 
not particularly significant for individual objects, but we dis- 
cuss the issue further because of the intriguing behavior of 
Avis and AvL ya with respect to one another and because of 
the greater significance of the result in higher S/N composite 
spectra discussed below. Three of the 1 1 galaxies with mea- 
sured Avis > also have Lya emission, and have (AvL yce ) = 
+708 ± 50 km s" 1 , ~ 250 km s" 1 higher than the average 
of the full sample; however, the average (AvL yQ -Avis) = 
622 ± 40 km s" 1 for this set of objects is nearly identical to 
that of the full sample. The relative consistency of the dif- 
ference AvL ya _ Avis in the Ha sample, as well as in much 
larger samples without the benefit of H a spectroscopy (e.g., 
Shapley etaD|200l LSteidel et al.ll200l . suggests a situation 
in which AvL ya moves in concert with Avis irrespective of 
whether the centroid of the absorption line velocities are blue- 
shifted with respect to systemic. 

Interestingly, 3 other galaxies of the 1 1 (see Fig. |2]i with 
Avis > are among those with spatially resolved velocity 
shear in the Ha emission line, meaning that the Ha spectrum 
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exhibits a si gnificant velocity offset as a function of spatial 
position; see lErb et alj (|2006c) for deta ils. Given the ov erall 
detection rate of shear in the sample of lErb et all (l2006cl) . we 
would expect to find ~ 1 such object in a sample of 1 1 . Fur- 
the rmore, of the 14 objects with tilted Ha emission lines in 
the lErb et al.ll200 6c sample, eight have high quality absorp- 
tion redshifts with (Avis) = -47 ± 35 km s _I ; 6 of 8 have 
Av IS > -60 km s" 1 . 

The apparent connection between decreased | Avis | and ob- 
servation of measurable velocity shear in Ha can be inter- 
preted in at least two ways. If the velocity shear indicates 
an unresolved merger, infalling gas (or an unusual amount 
of gas near zero velocity) could significantly reduce the av- 
erage blue-shift of the interstellar lines. The interstellar and 
nebular redshifts could also arise from differen t pieces of a 
merger , or even different galaxies altogether (see Quider et al. 
2009b for a possible example). Alternatively, if the tilted 
emission lines are caused by rotation that is most easily de- 
tected in nearly edge-on configurations, the outflows in these 
objects may be loosely collimated perpendicular to the disk 
as for many local starburst galaxies. If this were the case, 
one might expect to observe lower outflow velocities and per- 
haps stronger IS absorption near the galaxy systemic redshift 
for objects with smaller inclination angles to the line of sight. 
Both of these effects appear to play a role in the Na D IS line 
kinemati cs observed in a l arge sample of nearby star-forming 
galaxies (IChen et al.ll2.Q10t) . While possible inclination effects 
would not naturally account for highly redshifted Lya emis- 
sion, none of the objects with velocity shear identified from 
their Ha spectra happen to have Lya in emission. Inclination 
effects could be present in our sample, although the signifi- 
cant bulk outflows observed in the vast majority of the sample 
argue against collimation and projection effects being a major 
factor in most cases. 

As we will show in $4] below, we favor an explanation for 
many of the observed trends discussed in this section that 
hinges on the quantity of gas at or near zero velocity (and 
not on the overall outflow speed). 

3. THE RELATION BETWEEN BULK OUTFLOWS AND OTHER 
GALAXY PROPERTIES 

One of the advantages of the sample of z ~ 2 galaxies dis- 
cussed in this paper is that a large number of other galaxy 
properties are available to look for trends with respect to the 
bulk outflow properties. Most of the me asurements and i n- 
ferred quantities used here are tabulated in lErb et afll2006c llbl 
Among the parameters available are the Ha line widths <r v , 
the star formation rates inferred from the Ha line fluxes (cor- 
rected for extin ction according to the method outlined in 
lErb et al.ll2006bl) . the surface density of star formation Esfr, 
the dynamical mass Md y „ measured from a combination of a v 
and the observed physical size of the Ha emitting region, and 
the stellar mass M* inferred from SED fitting from the rest- 
UV to the rest optical/IR. The cold gas mass M gas is estimated 
by using the measured Ha surface brightness and galaxy size 
and a ssuming that the local Schmidt-Kennicutt dKennicuttl 
1998) relation between gas surface density and star formation 
rate applies. Finally, the total baryonic mass Mb ar = +M gas 
and the fraction of the inferred baryonic mass in the form of 
cold gas, u = M gas / Mbar, have been utilized. 

Table m summarizes the results of Spearman correlation 
tests between Avis, Avl yq , AvL yQ -Avis, and these other 
physical quantities. The number of galaxies available in the 
sample for each test, which is also given in table Q] varies 



TABLE 1 

Correlations of Av is and Av Lya with Galaxy 
Properties" 



Quantity 


-Avis 


Av Lya 


Av LyQ -Avis 




-2.08 (65) 


+1.81 (29) 


+0.74 (29) 


SFR C 


-1.52 (87) 


-0.04 (42) 


-0.08 (39) 


SsFR d 


+0.95 (81) 


+0.82 (37) 


+0.49 (35) 




-2.24 (57) 


+1.02(24) 


+0.14(24) 


M eas f 


-1.68 (73) 


+0.71 (36) 


+0.85 (34) 


M»s 


-1.93 (73) 


-0.37 (36) 


-0.48 (34) 


M bar h 


-2.66 (73) 


-0.10(36) 


-0.10(34) 




+1.72 (73) 


+0.64 (36) 


+0.93 (34) 



a All values are the number of standard deviations from the null 
hypothesis that the quantities are uncorrected, based on a Spear- 
man rank correlation test. Negative values indicate anti-correlations 
between the quantities. The number in parentheses following each 
value is the number of galaxies in the sample used to evaluate the 
correlation. 

b Velocity dispersion measured from the Ha emission line. 

c Star formation rate, in Mq yr , measured from the intensit y of the 

Ha emission line, and corrected for extinction as in Erb et al. 2006b. 

d Average star formation surface density, a s in Erb et al. 2006b]. 

c Dynamical mass, as tabulated in Erb et al. 2006a. 

f Cold gas mass, estimated fro m the star formatio n surface density 

and the observed Ha size, as in Erb et al. 2006a fl 

E Stellar mass, estimated from SED fitting, from Erb et al. 2006c 

h Total baryonic mass, M* +M gas . 

1 Gas fraction, M gas /M(, iu , as in Erb et al. 2006a c. 

depending on the quantity being evaluated. The tests have 
been conducted against the absolute value of the quantities 
Avis and AvLy Q so that the sense of any correlations is posi- 
tive when the bulk velocity differences increase with the other 
physical characteristic. None of the quantities considered is 
significantly correlated with Av LyQ or Av LyQ - Avis at more 
than the 95% (2a) confidence level, although the measured 
Ha velocity dispersion a v may be marginally correlated, in 
the sense that Avl yq is larger for objects with larger a v . 

The measured values of Avis, on the other hand, show 
greater than 2cr deviations from the null hypothesis (no cor- 
relation) for the quantities er v , Mdyn, and Mf, ar - While a v and 
Mdyn are correlated with one another by virt ue of the fact tha t 
the former is used to calculate the latter (see lErb et al.l l2006c). 
the two mass estimates Md y „ and Mt ar are measured using 
completely independent methods. lErb et al] d2006c) show 
that these mass estimates track each other very well, inspir- 
ing some confidence that both are providing reasonable esti- 
mates of the true masses within the Ha emitting regions of the 
galaxies. Within the Ha sample considered here, the inferred 
value of Mb ar ranges from ~ 1 x 10 10 to ~ 3 x 10 11 M , and 
was shown to track Mdyn much more tightly than , the stel- 
lar mass. Fig. [3] shows plots of Md yn and Mh ar versus Avl yq 
and A vis. The sense of the possible correlations of Avis with 
galaxy mass is that Avis (the centroid velocity of the inter- 
stellar absorption lines) tends to be closer to v = for objects 
with higher mass. The more significant of the mass corre- 
lations, with Mbar (significant at the 99% confidence level), 
appears to arise because none of the galaxies with baryonic 
masses smaller than ~ 3 x 10 10 M Q have measures of Avis 
near zero, while at masses Mbar > 3 x 10 10 Mq, a significant 
fraction do. 

There are no significant correlations between | Avis | and pa- 
rameters related to star formation. The lack of correlation 
is interesting in light of recent studies relating outflow ve- 
locity to galaxy properties in local starbursts; for example, 



The Circum-Galactic Medium at z = 2 — 3 



7 




-600 -400 -200 



ZOO 400 600 800 1000 
Av,„„ 



FIG. 3. — Plots showing the dependence of bulk outflow velocities on dy- 
namical mass estimated from Ha line widths and observed sizes (top panels) 
and total baryonic mass (i.e., the sum of inferred stellar and cold gas masses) 
(bottom panels). 



both lMartinl (120051) and lRupke eTaTl (120051) find evidence for 
correlations of outflow speed with SFR and dynamical mass 
(represented by the circular velocity), in the sense that the 
outflow velocities increase with both quantities (i.e., appar- 
ently the opposite to what is observed in the z — 2 sample). 
These relatively local galaxy samples have a dynamic range 
spanning ~ 4 orders of magnitude in SFR, and include dwarf 
starbursts with v c ~ 30 km s" 1 and SFR < 1 M© yr" 1 , i.e. a 
much larger range than present in our z ~ 2.3 sample. The 
observed correlations in the local sample flatten for galaxies 
with SFR > 10-100 M Q yr -1 dRupke et al.ll2005l) . the approx- 
imate range of SFRs in our current sample. In other words, 
no trends are present in the low-redshift samples when only 
galaxies with parameters characteristic of our z ~ 2 galaxies 
are considered. Also, as we discuss below, other effects on 
the kinematics measured from centroid velocities would likely 
mask the presence of a weak corre l ation i f it were present. 

At higher redshift, IWeiner et ail (120091) also found that the 
inferred outflow velocity in composite spectra of z ~ 1.4 
star-forming galaxies is a slowly increasing function of SFR, 
v om oc SFR 03 . However, because thei r composite s pectra 
are of relatively high spectral resolution. IWeiner et al.l (2009) 
measure Avis as the velocity at which the IS absorption 
reaches 90% of the continuum value, i.e. close to the max- 
imum blue shifted velocity rather than the centroid. They de- 
compose the Mg II IS line into a "symmetric" and "outflow- 
ing" component, and find that the strength of the symmetric 
component is very steeply dependent on stellar mass- in fact, 
it is much steeper than the variation of the outflowing compo- 
nent. As in our sample, the lowest-mass galaxies are consis- 
tent with having zero symmetric component (see their Fig 13 
and Table 1). 

Clearly, the centroid velocity of strongly saturated IS ab- 
sorption lines is in many ways a blunt tool for characterizing 
the velocity of outflowing gas. First, the lines may include 
a significant amount of interstellar gas at or near zero veloc- 
ity with respect to the galaxy systemic redshift, which would 
have the effect of decreasing the measured value of |Avis| 



even i f outflow velocitie s w ere substantial. W hile [Martin 
l2005t iRupke et al.l 120051 and IWeiner et al.l 120091 removed a 
component of interstellar absorption centered at zero veloc- 
ity before evaluating the velocity of outflowing material, it is 
generally not possible to do this for individual spectra in our 
z — 2 sample because of more limited spectral resolution. It 
is also not clear that subtracting a symmetric v = component 
of IS absorption is the best approach for evaluating the kine- 
matics of circumgalactic gas. It is possible, or perhaps even 
likely, that the IS line profiles would be sensitive to infalling 
material, which when seen in absorption against the galaxy's 
UV continuum would tend to be redshifted by up to a few hun- 
dred km s -1 with respect to a galaxy's systemic redshift; if so, 
the observed redshifted wing of IS absorption could be sub- 
stantially stronger than the blue-shifted port ion of the v ~ 
profile. As we have mentioned (see also jsl8.4l below). gas with 
positive (red-shifted) velocities with respect to z sys would be 
expected in at least a fraction of galaxies if infalling cool gas 
is present with a non-negligible covering fraction. 




-600 -400 -200 



200 400 600 800 1000 



FIG. 4. — Dependence of the Ha velocity dispersion ov on centroid veloci- 
ties AvLy a and Avrs measured from the far-UV spectra. 



Finally, we address whether outflows infl uence the Ha lin e 
widths used to calculate dynamical masses (lErb et al.ll2006cT) . 
In Fig. |4] we show a comparison of Avl yq and Avis with 
the Ha velocity dispersion <r v . As mentioned above, there 
is a marginally significant (~ 2a) tendency for |Avts| to be 
smaller in galaxies with larger values of a v . The correlation 
is likely to be related to the trends with galaxy mass we have 
already noted. If strong outflows were significantly influenc- 
ing the observed Ha lines, one might expect objects with the 
largest | Avis | also to have larger a v , a trend opposite to what 
is observed. As further evidence of the disassociation of the 
observable Ha emission and the winds, composite UV spectra 
(discussed in © constructed by combining individual spectra 
shifted to the nebular redshift show that Ha is at the same red- 
shift as the stars. Additionally, observations of local galaxies 
suggest that Ha emission from outflowing mate rial would fall 
far bel ow our detection threshold; for example. iLehnert et alj 
( 1999) studied the extended Ha emission from the superwind 
in the starburst galaxy M82, finding that it has a total lumi- 
nosity of 2.4 x 10 38 ergs s" 1 and comprises only ~ 0.3% of 
the total Ha flux. Since the typical S/N of our z ~ 2.3 Ha 
observations is ^10, it would be extremely difficult to rec- 
ognize a component of em ission coming from wind material. 
Finally. IColina et aTl (120051) have also found, through integral 
field spectroscopy of Ha emission in local ULIRGs, that the 
central velocity dispersions are unaffected by outflows. 

4. INFERENCES FROM COMPOSITE FAR-UV SPECTRA 
4.1. The use ofzna to measure galaxy systemic redshifts 
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The existence of a large sample of galaxies for which both 
rest-frame optical (Ha) and rest-frame far-UV spectra are 
available provides an opportunity to evaluate high S/N com- 
posite UV spectra formed from unusually precise knowledge 
of the systemic redshifts of the galaxies. A stacked compos- 
ite UV spectrum was formed by shifting each spectrum into 
the rest-frame using ZHa, scaling based on the flux density in 
the range 1300- 1500 A, re-sampling onto a common wave- 
length scale of 0.33 A per pixel, and averaging at each disper- 
sion point, with outlier rejection. Thus, the composite is an 
unweighted average of all 89 galaxies having both Ha spectra 
from NIRSPEC and high quality UV spectra from LRIS-B; 
the result is shown in Fig. [5] A measurement of the cen- 
troids of weak stellar photospheric absorption features (S V 
A1501.76, O IV A13 43.354, and C III Al 171.76 were used 
for this purpose; see IShaplev et alj|2003t [Pettini et al.ll2000l) 
in the rest-frame composite spectrum verifies that the Ha red- 
shifts are very close to the systemic redshifts of the stars in the 
galaxies, with a mean velocity of v* = 13 ± 24 km s -1 . Com- 
posites formed from various subsets of the data yield a simi- 
lar level of agreement: for example, a composite formed from 
the subset of 28 galaxies with significantly deeper LRIS-B 
spectra (with total integration times ranging from 5-20 hours 
compared to the typical 1 .5 hours) has v* = 2 ± 36 km s" 1 . 

4.2. Velocity widths of IS lines and the maximum outflow 
velocity 

One of the quantities that can be evaluated with increased 
confidence using the new sample is the velocity extent of the 
interstellar absorption lines and Lya emission, in addition to 
the centroid velocities discussed in the previous section. Of 
particular interest are 771s, the FWHM of the interstellar ab- 
sorption lines, and 77 LyQ , the FWHM of the Lya emission 
line, corrected for the instrumental resolution. Unfortunately, 
as discussed above, the instrumental resolution is not always 
known precisely, because it depends on the seeing-convolved 
size of the galaxy compared to the l'/2 slits used for all of the 
LRIS-B spectroscopy. For objects illuminating the slit uni- 
formly, the spectral resolution is measured to be 450 km s" 1 
(FWHM) for the 400 line/mm grism used for almost all of 
the optical spectra of galaxies in the Ha sample; however, for 
typical seeing of ~ 0.6-0.8" , galaxies in our sample have 
FWHM~ 1 .0" . Given this, we would expect that the actual 
spectral resolution is FWHM~ 370 km s _1 (R = 800). We have 
verified this for both the composite spectrum of all Ha galax- 
ies, and for a smaller subset of galaxies with very deep spec- 
troscopic integrations, by measuring the spatial size of each 
object in the slit direction. Assuming that the Ha redshift 
uncertainties are a z ~ 60 km s -1 and that individual spectra 
have cr res ~ 160 km s" 1 (i.e., FWHM/2.355), the estimated in- 
trinsic FWHM of IS absorption features and Lya emission 
in the composite spectrum are r]/s ~ 540 - 570 km s -1 and 
77L yct ~ 620-650 km s" 1 , respectively. 

We have recently obtained LRIS-B spectra of galaxies se- 
lected in the same way as the Ha sample, but observed using 
the 600 line/mm grism instead of the 400 line/mm grism used 
for the vast majority of the Ha sample. These observations 
provide spectral resolution 1 .68 times higher for the same slit 
width and object size. As a test of our ability to measure line 
profiles with marginally resolved data, we assembled a sam- 
ple of 102 spectra all obtained on the same observing run and 
with consistent observing conditions (seeing of ~ 0.6-0.7" 



FWHM evaluated at 4800 A) in a field which remained very 
close to the zenith at the Keck Observatory, thus minimizing 
issues of differential atmospheric refraction. We applied the 
rules given in equation 4 above to shift the spectra into the 
rest frame, and produced a composite spectrum, shown in the 
bottom panel of Fig. [5] with an effective spectral resolution 
of FWHM- 225 km s" 1 (R ~ 1330). The stellar photospheric 
absorption features in the stacked spectrum have v* = -2 ± 10 
km s" 1 , illustrating that the rules for estimating the systemic 
velocity from the UV spectra work very well on average, and 
that the higher resolution spectra are advantageous for pro- 
ducing more accurate wavelengths for weak (stellar and IS) 
features. 

The line profiles of the R = 800 and R = 1330 composite 
spectra are remarkably similar, for both Lya emission and the 
strong interstellar absorption features, as shown in Fig. [6] The 
spectra yield the same value of r]iy a and 771s after accounting 
for the difference in spectral resolution. One is particularly 
interested in the asymmetry of the line profiles and the value 
of the maximum blue-shifted velocity v max , which we define 
as the velocity at which the blue wing of the IS absorption 
lines meets the continuum. We find that v max , although diffi- 
cult to measure for typical spectra of individual galaxies due 
to limited S/N, is not strongly dependent on spectral resolu- 
tion; Figure [6] shows that \v max \ ~ 700-800 km s" 1 for both 
composite spectra. As discussed above, the C IV and (to a 
lesser extent) Si IV doublets in galaxy spectra have contribu- 
tions from both the IS lines and from the P-Cygni stellar wind 
lines from massive stars, and thus the IS component must be 
separated from the stellar feature in the process of fitting the 
local continuum. An example of a co ntinu um fit near the Si IV 
and C IV features is shown in Fig. 14.21 Fortunately, the P 
Cygni feature is generally both broader and shallower than 
the IS components of these lines, so that while the continuum 
uncertainties are larger than for unblended features, they do 
not prevent an accurate measurement in relatively high S/N 
composite spectra. 

The velocity profiles of some of the strongest spectral fea- 
tures in the stacked spectrum of thegalaxies observed with 
R = 1330 are over-plotted in Fig. [8] Evidently, \v max \ ~ 
700-800 km s _1 is a generic feature of the spectra of these 
rapidly star-forming galaxies, in spite of the fact that the av- 
erage centroid of the IS line profiles is more modest, with 
Avis ~ -165 km s -1 from the previous section. The apparent 
values of \v max \ are relatively insensitive to spectral resolution. 

The line profiles in the spectra of individual galaxies can, of 
course, vary considerably. Fig.|9]shows the spectra of two in- 
dividual galaxies to illustrate the point: one is MS1512-cB58, 
the z = 2.72 9 lensed LBG whose sp ectrum has been analyzed 
in detail bvlPettini et al.l {2000, 2002); the other is Q0000-D6 
(IShaplev et al.ll2003l) . a bright LBG at z = 2.966 observed at 
a comparable spectral resolution of ~ 1500. These spectra 
show clear differences in the details of the profiles and with 
the apparent covering fraction of the continuum, particularly 
for the low-ionization species which differ in apparent optical 
depth by a factor of ~ 2. Clearly, Lya emission is promi- 
nent in D6, but very weak in the spectrum of cB58 (see also 
lOuider et al.ll2009al) . The spectrum of Q0000-D6 also has an 
unusual high ionization component that produces clear Lyman 
a absorption in the apparent blue wing of the Lya emission 
line, as well as in the high ions (but is less prominent in the 
lower ions) near v = 0. Still, the velocity envelope for the blue- 
shifted material is remarkably consistent with that seen in the 
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FIG. 5. — Composite rest-frame far-UV spectra for two independent samples of z ~ 2.3 galaxies. The top panel is an average of the 89 spectra in the Ha 
sample, with R = 800, after normalizing each to the same relative intensity in the range 1300-1500 A. The bottom panel is a composite of 102 galaxy spectra 
obtained with higher spectral resolution (R = 1330), shifted into the rest frame using equations[2]and|4] and scaled as the first sample before averaging. 



composites presented above: maximum blue-shifted veloci- 
ties of \vmax ~ 800 km s" 1 , roughly independent of ionization 
level. 

4.3. Trends with Baryonic Mass 

Returning to the trends in the centroid velocities of the IS 
lines noted in the previous section, it is instructive to exam- 
ine the mean line profiles of composite spectra selected by 
implied baryonic mass Mbar, the parameter most significantly 
linked to the observed kinematics of the interstellar absorp- 
tion features. Figure [10] shows the comparison between the 
halves of the Ha sample that are above and below an inferred 
Mbar = 3.7 x 10 10 Mq, the sample median. There are clear 
differences in the Lya emission line strength, which is simi- 
lar to tha t seen f or sub-samples of different metallicity as in 
lErb et al.1 d2006al) . The peak of the Lya emission line profile 
is shifted by ~ +200 km s _1 (from ~ +400 to ~ +600 km s" 1 ) 
for the higher mass sub-sample relative to that of the lower- 
mass sub-sample. 14 The profiles of the low-ionization inter- 
stellar lines (Cll A1334 is in the cleanest spectral region and 

14 No significant correlation was found between Mbm and AvL yce in {fj] 
but many of the galaxies in the higher mass sub-sample had Lya emission 
that was too weak to measure, resulting in a very small sample. 



illustrates it best) may be indicating the root cause of the kine- 
matic trends discussed in |J3] the higher mass sample exhibits 
stronger IS absorption at or near v ~ 0, while the profiles are 
nearly identical in their behavior near v = -|v maf |. This "ex- 
cess" low velocity material in the higher-mass sub-sample - a 
shift of ~ 200 km s -1 in the red wing of the IS line profiles- 
systematically shifts the centroid of the IS velocity distribu- 
tion by ~ +100 km s -1 relative to the lower mass sub-sample, 
while the blue wing of the profile exhibits no clear trend with 
Mbar- The excess low-velocity material in the higher mass 
sub-sample is not obvious in the C IV absorption profile, for 
which the 2 profiles appear to be nearly identical for v 
(note that C IV was not used to measure Avis for any of the 
galaxies in the sample, because of the dependence of the rest- 
wavelength for the blend on the relative strength of the lines 
of the doublet). 

Fig [TT] shows the residual apparent optical depth At(v) for 
3 relatively isolated low-ionization transitions for the high- 
and low- M\, ar sub-samples. By this we mean the additional 
optical depth as a function of velocity that when added to the 
line profiles of the low-Mh ar sub-sample would produce line 
profiles identical to those of the high-M/^, sub-sample, i.e., 



4m(v)=4n(v)e 



-Ar(v) 



(5) 
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FIG. 6. — A comparison of the Lya and IS line profiles for the full Ha 
sample (red) and the composite of 102 identically-selected objects observed 
with higher spectral resolution (blue), where the latter have been shifted into 
the rest frame using the rules defined in i|2] F° r blended features such as 
O I A1302/Si II A 1304 and C IV AA1548, 1550, a separate panel is given for 
each; in both cases, one should look at the longer of the two components for 
the profile near v ~ and the shorter one to gauge the extent of the most 
blue-shifted velocities. Note the very similar velocity profiles for the two 
independent samples, for both the minimum and maximum velocities. 
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FIG. 7. — A portion of the composite spectrum shown in the lower panel of 
Fig.|5]in the vicinity of the Si IV and C IV doublets. The dashed red curve is 
an example continuum fit used to normalize the spectrum for measuring the 
strength of the IS component of these lines. Note that the continuum has been 
adjusted to remove the broad absorption due to the stellar wind component 
of C IV; a similar adjustment was made to remove the stellar component of 
Si IV, though it is much weaker than the C IV feature. The weak emission line 
within the broad C IV wind absorption feature is nebular Si II* A1533, one of 
several weak excit ed fine structure em ission lines observed in the spectra of 
z — 2-3 galaxies (Shapley et al. 2003). The weak absorption line near 1501 
A is photospheric S V. 

where //,,„ and //„, are the spectral intensity of the high-M^,, 
and low-Mt, a r sub-samples, respectively. The spectrum At(v) 



FIG . 8 . — The velocity profiles of strong IS lines and Ly a emission features 
relative to systemic for the composite UV spectrum of the R = 1330 sam- 
ple. The C IV profile is a solid blue histogram for the A 1548 component and 
dashed for A1550. 
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FIG. 9.— As for Fig.|6] comparing the spectra of two individual objects 
having particularly high quality UV spectra as well as accurate systemic red- 
shifts from stellar absorption features. The spectrum of MS1512-cB58 (the 
R ~ 1600 spectrum from Pettini et al. 2000) is in red, while the spectrum 
of Q0000-D6, obtained with R ~ 1300, is in blue. Despite significant dif- 
ferences in spectral "morphology" (e.g., D6 has strong Ly a line emission, 
while cB58 has strong absorption, and the ratio of the strengths of high ion- 
ization lines to low-ionization lines is quite different), the range of velocities 
spanned by the outflows are quite similar in both cases, with the maximum 
velocities of v ~ 800 km s~' . Note that the Sill A 1526 line of D6 is affected 
by absorption from another system in the blue wing of the profile. Both of 
these galaxies would be in the Mbar < 3.7 X 10 10 Mq (i.e, lower baryonic 
mass) sub-sample. 

has a peak at v ~ 0, a centroid at v = +154 km s" 1 , and a veloc- 
ity width of <r v ~ 120 km s" 1 after correcting for the effective 
instrumental resolution. The excess apparent optical depth 
Ar in these transitions accounts for ~ 25% of the equivalent 
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FIG. 10.— As for Fig.|6] comparing the composites of subsets of the Ha 
sample depending on their inferred baryonic mass M|, ar . The composite from 
objects with Mbar > 3.7 X 10 10 Mq is plotted in red, and that from objects 
with M(,ar < 3.7 X 10 10 Mq in blue. Note the differences in the profiles of the 
low-ions (e.g., Sill, OI, CII) for v 0, where the high baryonic mass objects 
appear to have stronger low-ion absorption at velocities closer to systemic. 
These differences do not seem to be present for the higher ionization lines 
(e.g., see the CIV A1550 panel), nor are there differences in the maximum 
blue-shifted velocities. 

width of the full low-ion profiles of the high-M^,,- sub-sample. 
The additional component of absorption in the high-M/, flr sub- 
sample also changes the average interstellar line width from 
771s = 540 km s -1 for the lower mass sub-sample to 771s = 660 
km s -1 for the higher mass sub-sample. 

We will return to a discussion of the possible implications 
of the baryonic mass dependence of the IS line kinematics in 



5. LYMAN a EMISSION 

5.1. Observed Trends 

We have shown in §|2] above that the centroid velocity of 
Lya emission AvL yQ exhibits larger scatter than Avis relative 
to the systemic redshift defined by Ha emission, but that the 
relative consistency of Avl Y q - Avis suggests a causal link 
between the two kinematic measures. In this section, we at- 
tempt to understand the nature of this relationship in greater 
detail. 

In U above, we presented evidence that the kinematic pro- 
files of strong interstellar lines in the spectra of rapidly star 
forming galaxies exhibit greater variation near v = than at 
large blue-shifted velocities. Moreover, the presence of sig- 
nificant low-ion absorption at v ^ strongly affects both the 
apparent velocity of Lya emission and the centroid velocity 
Avis of the IS lines, though v max , the maximum blue-shift, 
remains essentially unchanged. 

A commonly-adopted "toy model" (see e.g. iPettini et al.1 
2002; Adelberger et al. 2003) used to interpret the kinemat- 
ics of Lya emission and IS absorption lines in starburst 
galaxy spectra involves a roughly spherically-symmetric out- 
flow resulting in generally blue-shifted absorption as seen 
from earth. Because Lya photons resonantly scatter, they 
escape from the nebula only when they acquire a velocity 
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FIG. 1 1 . — The relative optical depth Thigh ~t\ow for selected low-ionization 
absorption lines, versus velocity relative to systemic, where these quantities 
refer to the composite spectra of the high-M[, al and low-M(, ar sub-samples. 
The histograms are for Sill A1260 (green), CII A1334 (red), and Sill A1526 
(magenta) interstellar line profiles, where the heavy black curve is the aver- 
age residual optical depth for the 3 transitions. Note the absence of signif- 
icant residuals at blue-shifted velocities (with the possible exception of lo- 
cal "peak" near v = -300, which represents slightly deeper absorption in the 
spectrum of the low-Mh iU sub-sample). The excess absorption in the Mh iU - 
high sub-sample has its centroid at v = +154 km s and its peak near v ~ 0. 



such that the optical depth to scattering in material which lies 
(physically) between the observer and the emitted Lya photon 
becomes small. Perhaps the easiest way for a Lya photon to 
reach earth is to acquire the velocity of outflowing material on 
the far side of the galaxy, but to be emitted in the observer's 
direction, so that the photon has been redshifted by several 
hundred km s" 1 relative to the bulk of the material through 
which it must pass to reach us. This picture would explain 
qualitatively why the dominant component of Lya emission 
always appears redshifted relative to the galaxy systemic ve- 
locity. Neglecting radiative transfer effects for the moment, 
if the range of velocities of outflowing gas is similar on the 
"far" side of the galaxy to what we observe (through the IS 
absorption lines) on the "near" side, then one would expect 
the maximum blue-shifted velocity \v max \ to be comparable to 
the maximum velocity observed in the red wing of the Lya 
emission line. One could then explain different spectral mor- 
phology near the Lya line by altering the distribution of Lya 
optical depth as a function of velocity for the material between 
the observer and the far side of the outflow. 

For example, in the simplest possible scenario, one could 
imagine that the difference in the Lya line profiles shown in 
Fig|T2](or any other set of Lya profiles) could be explained 
by altering the distribution of Tj_, ya (v) such that only the most- 
redshifted photons have an appreciable chance of making it 
through the intervening H I. It is then easy to see that increas- 
ing TLy a (v) near v = will shift the apparent peak and centroid 
of Lya emission toward the red, since Lya photons emitted 
from material with line-of-sight velocities near v = will have 
no chance to reach the observer- only those with the most ex- 
treme redshifted velocities will find low enough TL yQ to make 
it through in our direction. Figure [12] shows the Lya profiles 
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FIG. 12. — A comparis on o f the composite Lya emission profiles of the 
same samples as in Fig. |10| Also shown (green histogram) is the inferred 
function 1 + ArLyo, (v), where Ai\ ya (v) is the excess apparent optical depth 
in Lyct (as a function of velocity) required to produce the weaker of the two 
Lya profiles from the stronger. The black histogram reproduces the low-ion 
residuals from Fig. 1111 

for the two samples split by Mbar as in the previous section- 
the green histogram shows 1 + ArL yQ (v) with ArLy Q defined 
so as to produce the red profile (high Mb ar ) from the blue one 
(low M bm ), i.e.,. 



4m(v)=/; m (v)e 



-ATL yQ (v) 



(6) 



Here ArL yQ is the excess optical depth, over and above that 
already present for the low-Mb ar sub-sample; the effective op- 
tical depth, including both scattering and dust opacity, must be 
very significant even for the latter, since the equivalent width 
of the observed Lya emission line is ^20-40 times smaller 
than that expected for Case-B rec ombination and a n ormal 
population of high mass stars (e.g. JCharlot & Fall 199j]) . 

Fig. [12] also reproduces the excess optical depth 
At(C II,Si II) in the low-ionization metallic species from 
Fig. [TT|for comparison to Ar(Lya). While both have cen- 
troids with v > 0, Ar(Lya) peaks at v ~ +270 km s -1 , close 
to the velocity at which At(C II, Si II) begins to decrease going 
toward higher positive velocity. The differences are probably 
due to a combination of generally much higher optical depths 
in the Lya transition than for the low-ion metals (Lya pho- 
tons can only escape from regions having the lowest r(Lya)), 
as well as geometric effects relating to the possible location 
of gas seen in absorption (which must lie in front of the con- 
tinuum source) versus emission (which may lie either in front 
of or behind the continuum source). 

5.2. Understanding Lya Emission Line Kinematics 

The real situation is undoubtedly far more complicated. 
There have been several recent theoretical treatments of 
Lya radiative tra n sfer in t he con te xt of galactic outflows 
(IVerhamme et alJ 120061 120081 : ISchaerer & Verh amme 
20081: iHansen&Ohl 120061: iDii kstra et all l2006blal : 
Zheng & Miralda-Escude 2002), and even in the highly- 
idealized geometric configurations considered in these 
papers, a given line profile does not uniquely specify the 



combination of velocity field, H I column density, covering 
fraction, and dust opacity that applies to a given observation. 
For example, both bulk velocity fields and photon diffusion 
are capable of accounting for Lya photons that acquire 
large redshifts before escaping from a model galaxy. The 
profiles of interstellar absorption lines, particularly in the 
highest-quality spectra of individual objects, indicate that 
material with velocities ranging from ~ +200 to ~ -800 
km s" 1 exists in most observed ions for most galaxies (as 
discussed above). Essentially all of the ions observed in 
spectra of the quality presented here are strongly saturated, 
so that line profiles are best thought of as maps of covering 
fraction (hereinafter f c ) versus velocity. The simultaneous 
presence of neutral and singly ionized species along with 
higher-ionization species like C IV, with similar overall 
velocity envelopes, reinforces the idea that the ISM is a 
complex multi-phase medium. To make matters worse, we do 
not know where, in physical space, IS absorption at a given 
observed velocity actually arises. Since Lya photons must 
traverse this medium, experiencing typically thousands of 
scattering events, it is therefore extremely difficult to predict 
in detail what the emergen t Ly a profile will look like. 

IVerhamme et al.1 ((2006) and IVerhamme et alJ ((2008) con- 
sider a wide range of parameters, examining the effect on 
emergent Lya line profiles, for models of central (monochro- 
matic) point sources surrounded by an expanding shell of gas 
with varying H I column density and Doppler parameter (b). 
A generic prediction, as discussed above, is that Lya emis- 
sion (when present) will be very asymmetric, with the de- 
tails of the line shape depending on the assumed shell veloc- 
ity, Doppler parameter, and H I column density in the shell. 
For the expanding shell models, these authors predict that 
the peak of the Lya emission line should appear near ve- 
locity v ~ -2 x V a p, where V exp is the shell velocity, due to 
the combined effects of radiative transfer and the bulk veloc- 
ity of the scattering medium. The red wings of the Lya line 
are expected to be produced by photons scattering multiple 
times from the receding side of the expanding shell (as seen 
by an observer on earth); most of the migration of Lya pho- 
tons toward large redshifted velocities arises from absorption 
in the Lorentzian wings associated with high values of N(H I). 
Larger values of V exp and higher assumed values of N(HI) also 
accentuate the red wing of the Lya line and suppress regions 
close to the systemic redshift, moving both the peak and the 
centroid of Lya emission to higher velocities. A further pre- 
diction of the shell models is the presence of additional Lya 
peaks correspond to photons escaping in either the blue or the 
red wing of Lya associated with the approaching side of the 
shell; photons emitted from the red wing of the Lya line in 
the approaching side would then lie close to v = 0, while the 
blue wing would tend to form a peak with v —V exp . 

In a less-idealized situation, where the bulk velocity of the 
outflowing gas is not single-valued, but spans a more-or-less 
continuous velocity range of at least 800 km s" 1 , the dom- 
inant mechanism for the migration of Lya photons in both 
frequency space and real space becomes simpler, in some re- 
spects 15 . In a scenario with more gradual velocity gradients 
and a clumpy medium spread over a large range in galacto- 
centric distance r (rather than a "shell") associated with an 
outflow, photons can achieve velocities off-resonance by scat- 
tering their way through gas with a range of v; when they 

15 More general implications of this type of model are discussed in iUJ 
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escape the nebula and are observed as red-shifted Lya pho- 
tons by the observer, they would generally exhibit Doppler 
shifts that directly reflect the velocity with respect to systemic 
of the gas from which they were last scattered. Since the 
medium is clumpy, the column density N(H I) for individual 
"clumps" becomes less relevant than their velocity distribu- 
tion and covering fraction , since most Lya photons will scat- 
ter off the "surfaces" of clumps, rarely encountering regions 
where absorptio n or em i ssion in the broad wings of a line is 
important (see Neufeld (1991) for a discussion of this type 
situation). Note that the clumpy geometry is qualitatively dif- 
fer ent from, e.g., the "Hu bble expansion" model considered 
by Verhamme et al. ( 2006), which has v increasing smoothly 
with r. In this case, the radiative transfer is still important 
because the scattering medium is continuous and has associ- 
ated with it a particular N(H I) through which all Lya photons 
must pass on their way to larger r. The picture we are advo- 
cating has photons scattering from the r(Lya) ~ 1 surfaces 
of discrete clumps; this causes most Lya photons to emerge 
from scattering events near line center, and thus to acquire 
a Doppler shift characteristic of the velocity (with respect to 
systemic) of the most recent clump. In other words, bulk- 
motion induced velocity shifts, rather than radiative transfer 
effects, may be most responsible for the kinematics of the ob- 
served Lya emission line. 

In this context, the IS absorption lines provide a reasonable 
proxy for the velocity distribution of gas that will comprise 
the medium through which the Lya photons must scatter in 
order to escape in the observer's direction. 16 . If the flows are 
roughly spherically symmetric, it should be possible to seek 
consistency between the kinematics of Lya emission (which 
would probe the kinematics of gas on the receding side of the 
flow) and that of the IS absorption, which samples the blue- 
shifted, or approaching, side. The Lya photons will be most 
successful in escaping the galaxy when they acquire velocity 
shifts well off resonance of whatever material lies between the 
last scattering surface and the observer. This effect would tend 
to produce emission from both the redshifted and blue-shifted 
gas. 

Figure Qj] compares the velocity profiles of Lya emission 
and C II A1334 absorption for the same R = 800 and R = 1330 
composite spectra shown in Fig. [5] First, it does appear that 
both spectra exhibit a secondary blue peak in the Lya emis- 
sion line, with centroid velocity at v ~ -520 km s" 1 (although 
it is more evident in the spectrum with higher spectral resolu- 
tion), close to a "mirror image" (in velocity) of the dominant 
redshifted emission. While the strength of this secondary fea- 
ture in the composites is only ~ 5% of the primary redshifted 
component, it corresponds closely to the range of velocities 
seen in the most blue-shifted portion of the IS C II profile, 
with -350 v ^ -750 km s -1 . This range also corresponds 
to that over which the apparent optical depth of the IS absorp- 
tion is decreasing from its maximum, which extends in both 
composites over the velocity range -350 v km s -1 . 

5.3. A Simple Kinematic Model for IS Absorption and Lya 
emission 

In an attempt to produce simultaneously the salient features 
of the observed IS lines and Lya emission lines, we have 
constructed a very schematic kinematic model following the 

16 Note that the velocity field information provided by the IS lines has not 
been exploited in any of the models mentioned above. 
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FIG. 13. — Comparison of two independent galaxy sample composites il- 
lustrating the presence of a weak blue-shifted component of Lya centered at 
v = —520 km s -1 (indicated with a solid vertical line). The blue and green his- 
tograms are the velocity profiles of Lya and CII A 1334 for the R = 1330 
sample, while the red and magenta histograms are the same lines for the 
R = 800 Ha sample. The position of the secondary Lya peak corresponds 
to the maximum blue-shifted wing of the IS line profiles. 

line of reasoning outlined above, for the purposes of illustra- 
tion. In the models, we assume that optically thick gas (for ei- 
ther low-ionization metallic species or Lya) is present in two 
kinematic components, each of which is Gaussian in optical 
depth r for a given transition: one component is centered near 
v = and the other is outflowing with a velocity distribution 
(v) = -450 ± 150 km s" 1 . We choose 450 km s" 1 for the cen- 
troid of the outflowing component somewhat arbitrarily, but 
(in addition to producing line profiles that resemble the real 
ones) it is approximately equal to the escape velocity v esc at 
the virial radius of a dark matter halo of total mass ~9x 10 
Mq at z — 2.5, believ ed to be typical of the galaxies in our 
spectr oscopic sample dAdelberger et alj|2005bt IConrov et alj 
2008; see §8]). For the gas near the galaxy systemic redshift, 
we assume that the velocity distribution of the optically thick 
gas is similar to that of the H II regions trace d by Ha emission, 
i.e., er,,(Ha) - 100 km s" 1 (lErb et all2006bl) . All profiles have 
been convolved with an instrumental resolution of R = 1330. 

Fig- 03 shows example model line profiles for both IS ab- 
sorption and Lya emission, as seen by an observer on earth 
(note the righthand panel in each case is just a zoomed-out 
view). One needs to provide the normalization of the r dis- 
tribution (i.e., the maximum optical depth tq) for both com- 
ponents; in the first example, the v=0 component has a cen- 
tral optical depth of r ~ 20, and the outflowing component 
has a peak at r ~ 5 17 . The dark blue solid curve shows the 

17 The result is relatively insensitive to these numbers, so long as the lines 
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FIG. 14. — Examples of model spectra in which the velocity field of opti- 
cally thick gas in the galaxy is assumed to be comprised of two Gaussian (in 
optical depth r) components, both having a covering fraction of 0.6 (dashed 
profiles in each panel). The top 2 panels show the results for a model with 
one component outflowing with v'i = -450 ± 150 km s _1 , and an additional 
component with V2 = 0± 120 km s ; the bottom panels show a second ex- 
ample with i'2 = +150 ±200 km s -1 , with the velocity V] and optical depth 
of the outflowing component held fixed. Note how the kinematics of the ma- 
terial near v = influences the measured centroids of the net absorption lines 
(heavy blue curve). If the same gas kinematics are assumed for H I, the pre- 
dicted Lya emission profiles as seen by an observer would resemble the red 
curves; the Lya centroid is redshifted by an additional ^ 100 km s in the 
bottom model as compared to the top. . 

expected absorption line profile (assuming that the covering 
fraction f c = 0.6, for both components) for a saturated low- 
ion transition. The combination of the two assumed velocity 
components produces an absorption feature with a centroid 
(Avis) = -184 km s -1 and a rest equivalent width Wo ~ 2.2A- 
close to the typical values measured for the real galaxies. 
A second example, in the bottom panels of Fig. [14] has an 
identical outflowing component, but the v ~ component 
has been shifted to +150 km s" 1 (to approximately mimic 
the "excess" optical depth observed in the high-Mft ar sam- 
ple discussed above; see Fig. QT}, and broadened to a v = 200 
km s" 1 . The centroid of the resulting absorption feature is 
(Avis) = -22 km s" 1 , in spite of the fact that the outflowing 
component is identical by construction, with maximum blue- 
shifts of v max ~ -800 km s" 1 . 

The predicted Lya emission from the same galaxies is 
slightly more complex; for the moment we adopt a simple 
model in which the t(\v\) distribution is assumed to be the 
same as for the IS absorption lines. The probability of a scat- 
tering event in which a Lya photon is emitted by one atom 
and then absorbed by another with the same velocity (alterna- 
tively, one can think of it as the relative fraction of time a Lya 
photon spends resonantly trapped in gas with optical depth 
t(v)) is roughly e oc t 2 (v), since r(v) oc «h(v), where «h(v) is 
the density of neutral H atoms with velocity between v and 

are saturated 



v + dv.) The probability that a resonantly trapped Lya pho- 
ton escapes without being immediately scattered at or near 
the same velocity is oc exp[-r(v)]. Crudely, the probability 
that a Lya photon is scattered by an atom with velocity v and 
is able to escape in the direction of an observer without addi- 
tional scattering, as a function of v, is the Lya emission profile 

m, 



/(v)~Cr 2 (v) exp[-r(v)], 



(7) 



where C is a scale factor that can be adjusted to roughly ac- 
count for the destruction of Lya photons due to dust, which 
we assume is proportional to the total path length traversed 
by a photon as it works its way out of the nebula. Finally, 
we assume that the Lya photons reaching an observer who 
sees IS absorption profiles as in Fig. [14] have been scattered 
toward the observer with probability of "escape" according to 
eq.|7] To be clear, we do not concern ourselves with how Lya 
photons made it from the central H II regions to any part of the 
outflow- instead we simply assume that at any given moment, 
the relative probability for emission of a Lya photon is pro- 
portional to r 2 (v). We do not know the spatial distribution of 
the outflowing gas on the redshifted side of the flow as a func- 
tion of velocity, but we assume that it is spatially distinct (i.e., 
farther away from the observer) than the component of ISM 
near the galaxy systemic redshift, and that the blue-shifted 
component is closer. Further, we assume that the outflow ve- 
locity \v oM | is a monotonic function of galactocentric distance 
r on both the red- and blue-shifted sides of the flow. We found 
this assumption to be necessary to obtain kinematics and spec- 
tral morphology consistent with observations; it would also 
follow from the assumption that the line depth at a given ve- 
locity is directly related to / c , so that the highest velocity ma- 
teri al is farthest away fr o m the continuum source; see also §|7] 
and lWeiner eTaT] d2009h ; lMarfin & Bouchel ( 120091) . 

The red curves in the panels of Fig. [14] were produced in 
this way, and show the predicted Lya emission profile for a 
galaxy whose absorption profile (Lya or low-ionization metal 
lines) is given by the blue curve, as seen by the observer. It 
is worth noting that (to first order, in the absence of dust) the 
Lya emission profile is not strongly dependent on gas cov- 
ering fraction, since (as seen by an observer) the outflowing 
material is both the "source" and the "sink" for Lya photons, 
both of which would be reduced by the same geometric factor. 
We note a number of interesting features of the predicted pro- 
files: first, the centroids of Lya emission and IS absorption 
nicely reproduce those typically observed in the galaxy sam- 
ple; second, the extreme velocities of redshifted Lya emission 
and blue-shifted IS absorption reach 800 km s" 1 , as for the 
observed galaxies; third, the model in the top panel produces 
a blue-shifted component of Lya emission near v ~ -500 
km s" 1 - similar in both relative intensity and velocity relative 
to the redshifted component as in the composite spectra (see 
Fig.[T3Tl. In the model, the secondary feature arises from pho- 
tons scattered toward the observer from gas in the approach- 
ing part of the outflow; it is generally weaker than the red- 
shifted component if there is any dust in the galaxy because 
photons coming from the near-side of the outflow do not ap- 
pear as redshifted relative to the scattering medium because 
they are moving in the same direction. As a consequence, 
blue-shifted Lya photons have generally experienced a larger 
number of scatterings in order to escape in the direction of 
the observer, thereby increasing the probability of being ab- 
sorbed by dust. Finally, the apparent peak of Lya emission 
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is modulated primarily by the component of gas near v = 0- 
Lya emission is more redshifted, and weaker, when the ve- 
loci ty range spanned by the ISM in the galaxy is broader (see 
also lMas-Hesse et al.ll2003l) . 

In summary, the apparent velocity of Lya emission depends 
primarily on gas near the systemic redshift of the galaxy, 
while the extent of the red wing depends on the maximum 
velocity of material with appreciable column density and/or 
covering fraction f c . The extent of the red wing of Lya ap- 
pears to be rather consistent over sub-samples examined. The 
relative consistency of AvL yQ _ Avis reflects this dependency: 
as IS absorption lines shift toward more positive velocities, 
only the Lya photons that scatter off very high velocity mate- 
rial on the far side of the galaxy can penetrate the H I gas in 
the foreground. In any case, it appears that the spectral mor- 
phology of Lya and IS lines in rapidly star-forming galaxies 
may be more easily explained by the (observed)bulk veloci- 
ties of outflowing IS gas, as compared to more subtle radiative 
transfer effects. 

We return to a more physically-motivated outflow model in 

§13 

6. WHERE IS THE ABSORBING GAS? 
6.1. Mapping the Circum-Galactic Medium 

Lines of sight to star-forming galaxies can provide a great 
deal of information on the overall kinematics, chemical abun- 
dances, and (in some cases) estimates o f the mass flux of 
cool material entrained in a n outflow ("e.g. JPettini et alJuOOOl 
l2002t lOuider et al.1 l2009allbh . However, such observations 
contain little or no direct information about where the gas 
is located relative to the galaxy. A full appreciation for the 
physics of the outflows, and their effects on both the "host" 
galaxies and the local IGM, requires probing them along lines 
of sight from which the physical location of gas as a func- 
tion of ionization level and velocity can be disentangled. In 
practice, this is very difficult. 

In the absence of spatially-resolved spectroscopy, each 
galaxy spectrum provides only a surface-brightness weighted 
absorption profile in each observed transition, integrated 
along the entire line of sight from the "center" of the galaxy. 
Although ions of differing excitation appear to have similar 
velocity profiles when integrated along the line of sight, we 
do not have any way of knowing (e.g.) whether the bulk of 
the ClV ions at a particular v are distributed in the same gen- 
eral regions as those giving rise to C II absorption at the same 
velocity. The information available from the galaxy spec- 
tra themselves is necessarily crude given the relatively low 
S/N and spectral resolution of all but a handful of high red- 
shift galaxy spe ctra. The few high -quality spectra, of which 
MS1512-cB58 dPettini et alj|2002l) is perhaps the best exam- 
ple, show that the velocity profiles are complex, and that most 
of the lines accessible in lower quality spectra are strongly 
saturated. If so, then interpretation of crude measures of line 
strength such as equivalent width (or apparent optical depth 
vs. v) is ambiguous. The absorption arises in gas which may 
not be uniform across the entire continuum source. The ab- 
sorbing gas may be close to the galaxy, which is typically 
several kpc across, and some lines of sight through the (out- 
flowing) material may be completely absorbed, while oth- 
ers may be a l togeth er free of absorption. As discussed by 
Shapl evet alj ([2003), the equivalent width of a strong line is 
modulated primarily by the fraction of the continuum source 
covered by gas giving rise to the particular transition (which 



controls the depth of the lines relative to the continuum) and 
the velocity range in the foreground gas. 

Since the absorbing gas is almost certainly clumpy (e.g., 
the range of ionization level seen in absorption cannot coex- 
ist in the same gas, and yet they share a comparable overall 
velocity "envelope"), the apparent optical depth for saturated 
lines (as seen by an observer on Earth) arising from material 
at a given galactocentric radius r will depend on the charac- 
teristic physical scale of the "clumps" a c (r), and their num- 
ber density n c (r). The observed line strength (in the present 
case, the rest equivalent width Wo) will also depend on the 
distribution of line-of-sight velocity- large values of Wo re- 
quire a large line of sight velocity dispersion, the width of 
which depends on the sampled distribution of covering frac- 
tion. Clearly, any model for the outflows is under-constrained 
by the spect rum of the galaxy itself, even for high quality 
spectra (e.g- lPettini et al]l2002tlOuider et al.ll2009al) . 

However, high-resolution spectra of very closely-spaced 
sightlines toward gravitationally-lensed QSOs have provided 
constraints on the characteristic physical scales over which 
column densities of various transitions suc h as Lya, CIV, and 
lower-ionization species vary significantly (Rauc h et alJ 19991 
l200lUEllison et al.l l2004). These scales appear to be ~ 0.5-1 
kpc for C IV (thus, comparable to the "beam footprint" of the 
galaxy) and ~ 50- 100 pc for low ionization metallic species 
(much smaller than the typical beam footprint). If it is as- 
sumed that the QSO spectra are sampling metal-enriched gas 
similar to that found near galaxies like the ones in the cur- 
rent sample, these scales, together with measures of covering 
fraction, provide some useful information on the structure of 
the gas in the flows. Information on the distribution of gas 
transverse to the line of sight breaks some of the remaining 
degeneracies hindering a full phase-space (velocity, position) 
understanding of the circumgalactic gas. 

6.2. Galaxy-Galaxy Pair Samples 

One method for probing the geometry of outflows is to 
observe foreground galaxies found near the lines of sight to 
background QSOs, and to compare the kinematics of absorb- 
ing gas in the galaxy spectrum with that seen in the QSO 
spectrum. A picture of the structure of the gaseous enve- 
lope surrounding galaxies can then be built up statistically, 
by observing a sample of galaxies with a range of projected 
galactocentric distances from the background QSOs (e.g., 
lAdelberger et alj2005atlAdelberger et alj2003l) . This method 
has the advantage that one can choose the background QSOs 
to be very bright, so that very high S/N and spectral resolu- 
tion is possible, but the disadvantage is that such QSOs are 
rare, and one is "stuck" with whatever galaxies happen to lie 
within ~ 1 ' of the QSO. Of course, the smallest separations 
(e.g., 9 < 15" , or a projected physical galactocentric radius 
of 120 kpc ) will be especially rare. Ultimately, one is also 
limited by the quality of the (faint) galaxy spectrum if one is 
interested in a kinematic comparison using two different lines 
of sight through the galaxy. 

An alternative and complementary method is to use close 
angular pairs of galaxies at discrepant redshifts to vastly in- 
crease the sample sizes, at the expense of information quality 
(because both the foreground and the background objects are 
faint). An initial attempt to use this method was presented by 
lAdelberger et alj (12005a); here we investigate it using a sig- 
nificantly larger, higher quality, sample of galaxy pairs drawn 
from an ongoing d ensely-sampled spectroscopic survey (see 
Steidel et al. 2004 for early results; the full survey will be de- 
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scribed in detail elsewhere). 

The survey has been conducted over a total of 18 inde- 
pendent fields, 15 of which are centered on bright QSOs 
with zq — 2.6-2. 8 (h igh resolution spectra of the bright 
QSOs are used in jj6.5l >. The spectroscopic sample includes 
£ 2500 galaxy spectra in the redshift range 1.5 < z < 3.6, 
with (z) = 2.28 ±0.42. We used this galaxy survey catalog to 
identify pairs of spectroscopically-confirmed galaxies in each 
of 3 bins in angular separation: 8 < 5" , 5" < 9 <10" , and 
10" < 9 < 15" , which will be identified as samples PI, P2, 
and P3, respectively (see Table|2]i. At the mean redshift of the 
foreground galaxies of (z) = 2.2, the angular bins correspond 
to projected physical distances of b < 41 kpc, 41 < b < 83 
kpc, and 83 < b < 124 kpc for samples PI, P2, and P3, re- 
spectively (see Table |2). The conversion from angular sep- 
aration to physical impact parameter changes by only ±2% 
over the full redshift range of the foreground galaxy samples, 
so that we will use angular separation and impact parameter 
interchangeably. 

Within each pair sample, a number of additional criteria 
were imposed: 1) both galaxy spectra within the pair must 
have accurately-determined redshifts; the redshift-space sep- 
aration of the two galaxies in each pair must be > 3000 
km s" 1 to ensure that their geometry is unambiguous (i.e., so 
that we know which galaxy is "behind" the other); 2) the 
redshift separation Zb g -Zf g < 1.0 to ensure that the spectra 
contain a significant rest wavelength interval in common af- 
ter both are shifted to the rest frame of the foreground object. 
Redshifts were assigned to each galaxy using their Ha red- 
shift if available, and otherwise using equations 2 or 4, de- 
pending on spectral morphology. As discussed in $2.31 the 
uncertainty in the redshifts for the bulk of the galaxies in the 
samples are expected to be ~ 125 km s" 1 . After culling, the 
final pair samples include 42, 164, and 306 pairs for samples 
PI, P2, and P3, respectively. The properties of these three 
angular pair samples are summarized in Table [2] 

6.3. A Case Study: GWS-BX201 and GWS-BM115 

Among these pair samples, there is a handful for which the 
spectra were of sufficient quality to be analyzed individually. 
One exam ple, illustrated in Figs . [TBI and Fig. Q~6] (see also 
Fig. 4 of Adelberge r et aLll2005al) . allows us to compare the 
interstellar absorption lines in the spectrum of the z = 1.6065 
galaxy GWS-BM115 with the same features observed in the 
spectrum of a background galaxy projected only 1"9 away 
on the plane of the sky, GWS-BX201 (z = 2.173). In this 
case, the systemi c redshift of GW S-BM115 is known from 
Ha spectroscopy (lErb et alJl200 6b). and the optical spectra of 
both objects were obtained with spectral resolution R ~ 1500 
(i.e., ~ 2 times higher than most of the other spectra in the 
sample). These galaxies have apparent magnitudes g' = 23.5 
(BX201) and g' = 23.7 (BM1 15), ~ 1 magnitude brighter than 
typical galaxies in our sample. Using the systemic redshift of 
the foreground galaxy (GWS-BM1 15) measured from its Ha 
line, both spectra were shifted into the foreground galaxy's 
rest-frame, and continuum normalized, to produce the com- 
parison in Fig. [16] which shows the Al II A1670 and C IV 
AA1548,1550 transitions 18 . A summary of relevant measure- 
ments is given in Table [3] 

18 Because of the low redshift of BM1 15, there are fewer lines in common 
for these two galaxies than for most of the pairs in the 3 samples; the tran- 
sitions shown are representative of low-ionization and high-ionization lines, 
respectively. 




FIG. 15.— HST/ACS image of the GWS-BM1 15/BX201 galaxy pair, in the 
F814W filter. BM115 has z sys = 1.6065, based on the Ha line observed with 
Keck/NIRSPEC. BX201 lies 1'.'9 away on the plane of the sky, corres pond ing 
to an impact parameter at z = 1.6065 of ~ 16.1 kpc (physical). Fig. ll6l and 
Table [3] compare the interstellar lines observed in the spectrum of the fore- 
ground galaxy BM115 and those observed in the spectrum of BX201 at the 
redshift of BM115. 

There are several points worth making for this particular 
case: first, the lines seen in the spectrum of the foreground 
galaxy have blue-shifted centroids, with Avis ~ —192 km s -1 , 
which we have shown is typical of the objects in the Ha sam- 
ple. Although the C IV lines are very broad (FWHM~ 400 
km s" 1 after correcting for the instrumental resolution), one 
can still discern a separation of the two components of the 
doublet. The line of sight to the background galaxy (GWS- 
BX201) passes the foreground galaxy at a projected physi- 
cal distance (hereinafter "impact parameter", b) of b = 16.1 
kpc. Along this line of sight, the centroids of Al II and C IV 
are consistent with Avis = relative to the foreground galaxy 
GWS-BM1 15- fitting the two lines of the C IV doublet and the 
Al II feature yields a velocity offset Avis = ± 15 km s -1 and 
velocity width FWHM~ 300 km s" 1 after accounting for the 
instrumental resolution. For both transitions, the "offset" line 
of sight produces lines which are significantly stronger (in 
terms of equivalent width) than those seen directly toward 
BM115 itself (see Table [3}! Of course, the path length 
through the galaxy is different in the two cases; one line of 
sight passes through only the "front half" of whatever gas is 
associated with the foreground galaxy, while the sightline to 
the background object, although spatially offset, samples ma- 
terial both "behind" and "in front of" BM115 as seen by an 
observer on earth. The fact that the centroid of the absorption 
has Avis ~ is consistent with any picture in which material 
is distributed reasonably symmetrically around BM115; e.g., 
this would be expected for axisymmetric radial outflow from 
(or infall onto) BM115, reaching a galactocentric distance of 
at least r g > 16.1 kpc. 

As mentioned above, strong lines such as Al II A 1670 and 
C IV are very likely to be saturated; they clearly exhibit sat- 
uration along the b = 16 kpc line of sight, but the lines in the 
foreground galaxy spectrum itself do not reach zero intensity. 
Al II A 1670 and both components of the C IV doublet have 
approximately the same residual intensity of ~ 0.5 (Fig. IT6b. 
The inference of f c < 1 for cool outflowing material is not 
unusual, and in fact a covering fractio n of ~ 50% is quite 
typical of z ~ 2 - 3 galaxies (see, e.g., IShap lev et al.j 120031 ; 
lOuider et al.ll2009al) . 

Based on the HST/ACS image (I band) shown in FigurefTBl 
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TABLE 2 

Foreground/Background Galaxy Pair Statistics" 



Sample 


6 Range b 


b Range c 


Number 




(b) 


<%> d 




PI 


<5" 


10-41 kpc 


42 


3"8±0"8 


31 ± 7 kpc 


2.21 ±0.32 


2.64 ±0.36 


P2 


5" -10" 


41-82kpc 


164 


7'.'6±1'.'5 


63 ±12 kpc 


2.18±0.36 


2.64 ±0.41 


P3 


10" -15" 


82-125 kpc 


306 


12"5±1"5 


103 ±12 kpc 


2.14±0.33 


2.65 ±0.39 



il Galaxy pairs in the spectroscopic sample for which both the foreground and background galaxies have unam- 
biguously measured redshifts which differ by > 3000 km s _1 . 
b Range of angular separation of the galaxies, in arc seconds. 
c Range of impact parameter, in physical kpc. 

d Mean redshift for the foreground galaxies in pairs with the specified range of 8. 
c Mean redshift for the background galaxies in pairs with the specified range of 9. 



TABLE 3 

Absorption Lines at z fg in GWS-BX201/BM1 15 Pair 3 



Spectrum (b) W (C IV 1549) b W a (Al II 1670) Av IS c FWHM d 

(kpc) (A) (A) kms-' km s~' 

GWS-BM115 (z= 1.6065) 0.0 2.78 1.11 -192 ± 15 400 

GWS-BX201 (@z = 1.6065) 16.1 3.53 2.20 0± 18 270 



a For this galaxy pair, the spectra were both shifted to the rest frame of the foreground galaxy (BM 1 1 5) 
and then evaluated with respect to z = 1.6065. 

b Values include both components of the ClV doublet; the doublet ratios are 1.20 and 1.05, respec- 
tively. 

c Velocity of IS line centroids relative to z = 1.6065. 

d Velocity widths of IS lines, corrected for instrumental resolution 
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FIG. 16. — Comparison of the kinematics of representative interstellar ab- 
sorption lines in the spectrum of the z = 1.6065 galaxy BM115 (red) to the 
absorption in the spectrum of the background galaxy BX201 (z = 2.173, 1'.'9 
on the plane of the sky, or a projected transverse separation of ~ 16 kpc at 
z = 1.607), shifted into the rest frame of BM115 using the same systemic 
redshift (blue). The rest wavelengths of the transitions shown are indicated 
(the C IV profiles are plotted with respect to v = for the A1548.20 line). 
Note that in both cases the absorption in the background galaxy spectrum is 
stronger and shifted by ~ +190 km s compared to the features in the fore- 
ground galaxy spectrum, placing the centroid very close to the foreground 
galaxy systemic redshift. 



although the 2 galaxies have similar apparent magnitudes they 
have quite different surface brightness distributions: BX201 
has FWHM ~ 0.22" (~ 1.8 kpc, physical), while the UV 
light of BM115 has FWHM- 0.40" (~ 3.4 kpc, physical). 
Thus, the spectrum of BM115 itself indicates that, averaged 
over an area of ~ 10 kpc 2 centered at b = (assuming that the 
intensity-weighted diameter of the light distribution is given 
approximately by the FWHM), but through only "half" of the 
outflow, only ~ 50% of the "beam" area intercepts apprecia- 
ble outflowing low and intermediate ionization gas. The "off- 
set" line of sight (b = 16.1 kpc) produces absorption reaching 
zero intensity in both relatively low (Al II) and relatively high 
(C IV) ionization gas, indicating f c ~ 1 when averaged on a 
scale of ~ 2.5 kpc 2 , at a galactocentric radius of ~ 16.1 kpc. 

It is difficult to know from observations of a single galaxy 
pair whether differences in the kinematics and/or depth of ab- 
sorption lines arise from averaging over a smaller spatial re- 
gion (e.g., a more compact background source might subtend 
a region within which the gas is optically thick everywhere) or 
from a changing gas-phase covering fraction with impact pa- 
rameter b (such that the coherence scale of absorbing gas may 
be physically larger or smaller). This degeneracy is one of the 
disadvantages of resolved background sources. However, as 
we now demonstrate, it is possible to use a large ensemble of 
pairs to construct a physical picture of the size and structure 
of an average galaxy's circumgalactic gas distribution, again 
using stacked composite spectra. 

6.4. Composite Spectra of Galaxy-Galaxy Pairs 

Unfortunately, most of the galaxy pairs in the current spec- 
troscopic sample are considerably fainter than the GWS- 
BX201/BM1 15 pair, and were observed with somewhat lower 
spectral resolution. Nevertheless, the large size of the pair 
samples, coupled with the improved accuracy of systemic 
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FIG. 17. — Comparison of the absorption lines in composite spectra formed from galaxy-galaxy pairs. The top spectrum is the composite of all foreground 
galaxies in pairs from sample P3 with redshift differences > 3000 km s . The second spectrum is the composite of all background galaxies in pairs with 
angular separation A9 < 5" [sample Pl(bg)], formed by shifting each to the redshift of the foreground galaxy and averaging, as described in the text. The third 
and fourth spectra were formed the same way, for the background galaxies in pairs with 5-10" and 10-15" angular separations [samples P2(bg) and P3(bg)], 
respectively. 



redshifts estimated from far-UV spectral features (most of 
the galaxies in pairs do not have Ha measurements), allow 
for interesting measurements from composite spectra. For 
each angular pair of galaxies, the spectrum of the foreground 
galaxy in each pair was shifted to its own rest frame using 
the measured values of Avis and the rules given in eqs. 2 
or 4, and continuum-normalized. The resulting foreground 
spectra were then averaged within each pair sample, produc- 
ing stacked spectra which we refer to as Pl(fg), P2(fg), and 
P3(fg), respectively. A similar approach was used to produce 
stacked spectra of the background galaxies: each individual 
background galaxy spectrum was shifted into the foreground 
galaxy's rest frame using the same systemic redshift applied 
to the foreground galaxy spectrum, and continuum normal- 
ized. The strong IS absorption features at the redshift of each 
background galaxy constitute a source of noise that can af- 
fect the composite spectrum of the background galaxies (after 
shifting the spectra to z/ s , the lines will appear at different 
wavelengths for each spectrum). To minimize this effect, the 
strongest interstellar lines in each background galaxy spec- 
trum (at the galaxy redshift, Zbg) were masked when produc- 
ing the average spectrum at z/ g - The resulting composites 
Pl(bg), P2(bg), and P3(bg) are shown in Figure [17] together 



with (for display purposes only; in the analysis below we use a 
distinct stack of the foreground spectra within each pair sam- 
ple) the spectrum of the average of all foreground galaxy spec- 
tra in sample P3. 

Figures [18] [19] and [20] show line profile comparisons of 
Pl(bg) and Pl(fg), P2(bg) and P2(fg), and P3(bg) and P3(fg), 
respectively, for selected transitions. Table [4] summarizes 
the measured rest-frame line equivalent widths for the same 
stacked spectra. The quoted errors on the equivalent widths 
in Table [4] account for sample variance, continuum uncertain- 
ties, and measurement errors based on repeated measurements 
with spectra formed from subsets of the data. 

As might be expected, the composites Pl(fg), P2(fg), and 
P3(fg) are very similar, since they are averages of galaxy 
spectra for 3 sets of galaxies selected using identical criteria, 
with very similar redshift ranges (see Tabled. Nevertheless, 
they are completely independent, having been formed from 
the spectra of distinct samples of galaxies. The differences in 
measured line strength for the foreground galaxy composites 
indicate the approximate level of sample variance (< 10% for 
the strongest absorption lines). It is of interest to predict what 
the absorption line profiles of the spectra of the foreground 
galaxies would look like if the line of sight were to probe 
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TABLE 4 

Absorption Line Strengths at z fg in Galaxy Pairs" 



Sample 


(b) 


Lya 


Si II 1260 


CII 1334 


Si IV 1393 


Sill 1526 


C IV 1549 b 


Al II 1670 


Pl(fg) 







1.38 ±0.12 


1.79±0.15 


1.20±0.12 


1.34±0.12 


1.95 ±0.15 


0.97 ±0.10 




0° 


(4.7±1.0) d 


(2.01 ±0.18) 


(2.61 ±0.22) 


(2.04 ± 0.20) 


(1.96±0.18) 


(3.90 ±0.30) 


(1.42 ±0.15) 


Pl(bg) 


31 


2.01 ±0.15 


0.42 ±0.06 


0.90 ±0.08 


0.39 ±0.08 


0.37 ±0.06 


2.13 ±0.15 


0.40 ±0.08 


P2(fg) 







1.42 ±0.10 


1.74±0.15 


1.12±0.12 


1.52±0.12 


1.90±0.10 


1.10±0.12 




o c 


(4.8±1.0) d 


(2.07 ±0.15) 


(2.54 ±0.22) 


(1.90 ±0.20) 


(2.22±0.18) 


(3.80 ±0.20) 


(1.61 ±0.18) 


P2(bg) 


63 


1.23 ±0.20 


0.41 ±0.09 


0.67±0.12 


0.19±0.08 


<0.15 


1.18 ±0.15 


<0.20 


P3(fg) 







1.40 ±0.10 


1.62±0.10 


1.10±0.12 


1.41 ±0.10 


1.90±0.10 


1.01 ±0.11 




C 


(4.9±1.0) d 


(2.04 ±0.15) 


(2.37±0.15) 


(1.87 ±0.20) 


(2.06 ±0.18) 


(3.80 ±0.20) 


(1.47 ±0.22) 


P3(bg) 


103 


0.92 ±0.12 


<0.05 


<0.12 


0.12±0.06 


<0.04 


0.13 ±0.05 


<0.10 



a For each pair sample, the rest-frame equivalent width (in A ) for the stack of the foreground galaxy spectra are given in the first 
row; the results from the composite of the background spectra, shifted to rest frame of the foreground galaxy, are in the second row. 
b Values include both components of the C IV doublet. 

c Rest equivalent widths after applying corrections (see text) to represent sighttines through the entire galaxy at b = 
d Values for Lycv at b = are estimated from the observed strength of Ly/3 absorption. 
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FIG. 18. — A comparison of the absorption profiles in the composite spectra 
of sample Pl(fg) [red] and Pl(bg) [blue]. For each of the 42 parrs, the same 
foreground galaxy systemic redshift was used to shift both the foreground 
and background galaxy spectra to Zf g . Thus, the red curve represents the 
average galaxy absorption line spectrum (i.e., b = 0) while the blue curve 
is the average spectrum of the same ensemble ofgalaxies at mean impact 
parameter of (b) = 31 kpc (see Tableland Table [4] for sample descriptions 
and statistics. Note that the rest wavelength of the C IV doublet blend has 
assumed W (1548)/Wb(1550) = 1.4.) 

the full galaxy rather than only the part in the foreground. If 
one assumes that the strong lines are saturated, and that the 
kinematics of outflowing gas are similar on the far side of 
the galaxy as observed on the near side, the full line profile 
should include a reflection of the portion of the line profile 
with v < about v = 0. Generally, accounting for this effect 
increases the equivalent widths by a factor of 1.45 for low 
ionization species (e.g., CII, Sill), 1.70 for Si IV, and ~ 2 for 
C IV. These factors vary because of the differing strength of 
the absorption near v = - typically low-ionization species 
have stronger v ~ absorption, while for C IV there is very 
little. The corrected values of the line strengths for the fore- 
ground galaxy spectra are also indicated in Table |4j in the 
discussion below, we adopt the corrected values for the b ~ 



FIG. 19. — Same as for Figure [T8l for sample P2(fg) (red) and P2(bg) 
([blue). In this case, the blue spectrum represents the average absorption 
profile at impact parameters b = 4 1 — 82 kpc, for the same 1 62 galaxies which 
comprise the mean galaxy spectrum in red. (see Tableland table|4]for sam- 
ple descriptions and statistics.) 

sightlines. 

The comparison in Fig.[18]thus shows the difference in ab- 
sorption line profiles for lines of sight at b = 0, and those offset 
by V.5-5" ({9} = 3'.'8, or (b) ~ 30 kpc), for precisely the same 
set of 42 galaxies. As in the example of GWS-BM1 15/BX201 
above, the absorption line centroids in the (b) ~ 30 kpc line 
of sight are close to the foreground galaxy systemic redshift. 
For instance, the observed Lya absorption line in the Pl(bg) 
spectrum has a measured wavelength of 1215.62 A, only 12 
km s" 1 from v = 0; the estimated error in the line centroid from 
propagation of the uncertainties in z sys for individual objects 
is er(Av) w 20-25 km s" 1 . The other lines in the Pl(bg) spec- 
trum show similar behavior, in the sense that their centroids 
are much closer to v = than the same lines in the Pl(fg) spec- 
trum, which have v ~ -150 to -200 km s" 1 . It is also clear 
from Tableland Figure [l8]that the strength of low and inter- 
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FIG. 20.— Same as Figure[l8] for sample P3(fg) (red) and P3(bg) (blue). 
In this case, the blue spectrum represents the average absorption profile at 
impact parameters b = 82 — 125 kpc, for the same 306 galaxies which com- 
prise the mean galaxy spectrum in red. (see Table [2] and table [4] for sample 
descriptions and statistics.) 

mediate ionization absorption lines is significantly weaker for 
lines of sight with an average impact parameter of (b) = 3 1 
kpc (physical). Specifically, the strength of C II A1334 is re- 
duced by a factor of ~ 3, and SilV A1393 by more than a factor 
of 5, for sample Pl(bg) ((b) = 31 kpc) relative to Pl(fg). In- 
terestingly, the strength of the C IV absorption falls much less 
steeply in Pl(bg) versus Pl(fg), declining by only a factor of 
-1.8. 

The same general trend continues in the comparison of 
P2(fg) and P2(bg) (Figure [H and Table U, where P2(bg) 
represents an average spectrum for impact parameters b = 
41 -82 kpc with (b) = 63 kpc. P3(bg), sensitive to absorp- 
tion from gas with b = 82- 125 kpc ((b) = 103 kpc), has 
reached the point where the composite spectra formed from 
the low-resolution, low S/N galaxy spectra are just barely 
sensitive enough to detect absorption lines other than Lya . 
C IV absorption is detected with only marginal significance, 
and the low-ionization metal lines fall below the threshold of 
detectability. Interestingly, there is no significant change in 
equivalent width for the Lya absorption line between sample 
P2 and P3: both have Wo(Lya) reduced by only a factor of 
~ 2 relative to sample Pl(bg). We will discuss possible ex- 
planations for this behavior in §|7]below. 

A graphical summary of the contents of Table [4] is shown 
in Fig. |2U Points representing lines of sight with b ~ [i.e., 
Pl(fg), P2(fg), and P3(fg)] have been averaged, the result of 
which is plotted near b = 0. The small error bars show that 
the scatter among the independent sub-samples of foreground 
galaxy spectra is small. Ly/3 has been used in lieu of Lya for 
a measure of the strength of H I absorption in the foreground 
galaxy spectra, because of the contamination of Lya absorp- 
tion by Lya emission. Because both lines probably fall on the 
flat part of the curve of growth, the strength of Lya absorption 
is likely to be comparable to that of Ly/3. 

In spite of the limited sensitivity to weak absorption lines 
in the composite spectra, there are several trends worth not- 
ing. The strength of Si IV and both Si II lines track each 



other very closely out to b ~ 40 kpc. This is noteworthy be- 
cause it is to be expected when all 3 lines are strongly satu- 
rated and sample the same range of velocities; in particular, 
the two measured lines of Sill (A1260, A1526) have oscil- 
lator strengths / = 1.18 and / = 0.133, respectively. For a 
given column density N(Sill) the ratio of A 2 / for the two lines 
would differ by a factor of ~ 6.1; the fact that the equivalent 
widths of these two lines are equal out to ~ 40 kpc indeed 
supports strong saturation. This suggests the trend toward 
weaker absorption with increasing b is dominated by either 
a decreasing velocity spread or a decreasing f c ; the galaxy 
spectra are not of high enough resolution to distinguish be- 
tween the two. By (b) ~ 63 kpc, however, the measured Sill 
A 1526 line strength has decreased significantly with respect 
to Sill A 1260 (it is at least 3 times weaker) indicating that the 
A 1526 line is becoming unsaturated. The equivalent width ra- 
tio Wo(A1526)/Wo(A1260) implies an average total Si II col- 
umn density at (b) ~ 63 kpc of log N(Sill)~ 13.3 cm" 2 . In 
principle the two lines of the SilV doublet should allow an 
assessment of the degree of saturation, but the SilV A 1402 ab- 
sorption line is sufficiently weak in the composite spectra that 
it cannot be measured with useful precision. 

The other lines in Fig.|2T|exhibit trends similar to that of the 
Si lines, albeit with a shallower decline in line strength with 
increasing impact parameter. The equivalent widths for all 
of the ions in Fig. [21] fall precipitously somewhere between 
b ~ 60 and b ~ 90 kpc, with the exception of Lya. All of 
the species decline with galactocentric radius (or impact pa- 
rameter) as Wo oc b~P with 0.2 < j3 ^ 0.6 out to b ~ 50 kpc, 
whereas f3 £ 1 at larger b (again, Lya excepted). Because of 
line saturation, and the fact that only relatively strong lines 
can be detected using low-resolution spectra of faint galaxies, 
we know that the trends seen in Fig. yT|are dominated by a 
declining value of f c with galactocentric radius, along with 
a smaller range of line-of-sight velocities with increasing b. 
Schematic m odel s which reproduce the curves in Fig. 1211 are 
discussed in jj7.ll 



6.5. Lya Absorption at Larger Impact Parameter 

In order to trace the average Lya absorption line strength to 
larger galactocentric radii (b > 125 kpc), we made use of the 
S/N ~ 40- 100, high resolution (R ~ 30000) spectra of bright 
QSO s lying in 15 of the galaxy survey fields introduced in 
36.21 The QSO spectra were obtained using the HIRES spec- 
trograph on the Keck 1 telescope, with the UV cross-disperser 
providing simultaneous wavelength coverage from 3 100-6000 
A with very high throughput. For the time being, we degrade 
these QSO spectra (using only their high S/N) in order to ob- 
tain absorption line measurements analogous to those using 
the galaxy pairs above. As for the galaxy pairs, foreground 
galaxy redshifts were assigned based on equations 2 or 4, and 
the spectral regions of the QSO spectra corresponding to the 
Lya transition at the redshift of the foreground galaxy were 
extracted. The foreground galaxies were chosen to lie in two 
impact parameter bins: 125 < b < 200 kpc and 200 < b < 280 
kpc. Within each bin, 21 QSO spectral regions were shifted 
into the rest frame of 21 distinct foreground galaxies, i.e. so 
that v = lies at the galaxy redshift. The spectra within each 
bin were then averaged to obtain the mean Lya absorption 
profile. The resulting composites were then continuum nor- 
malized 19 and averaged to obtain the mean Lya absorption 

19 The strength of the absorption associated with the galaxies is measured 
relative to a new effective continuum level reduced by ~ 20%, the average 
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FIG. 21. — The dependence of IS absorption line strengths on the galactocentric impact parameter in physical kpc, for the galaxy-galaxy pair samples outlined 
in Tables [2] and [4] The values near b ~ are an average of the points for the 3 distinct foreground galaxy samples Pl(fg), P2(fg), and P3(fg). These values 
have each been corrected upward to account for the estimated contribution to the line strength from the "far side" of the gas distribution, assuming symmetric 
kinematics of outflowing material. Each color (for points and connecting line segments) represents a different ISM transition as indicated in the box legend; 
downward arrows on points indicate upper limits. The dashed curves using the same color coding are predictions for Wo(b) using the model described in the text 
(summarized in Table|5j. The two points for Lya represented by open squares were measured from HIRES spectra, as described in the text. The corresponding 
spectra are plotted in Fig. 1221 



profile. Although the background QSOs are point sources, 
and therefore not suitable for measuring gas covering frac- 
tion in a single line of sight, the ensemble of QSO sight- 
lines should provide an average absorption profile equivalent 
to what would be obtained using resolved background galax- 
ies. The two new points have been added to Fig. [21] to extend 
the measurement to b = 280 kpc, or an angular separation of 
~ 34" at (z) = 2.2. The averaged spectra in the galaxy rest 
frame for the two bins are shown in Fig. [22j the centroids of 
the resulting Lya absorption features are remarkably close to 
v = 0, another indication that th e estimated galaxy reds hifts 
are relatively free of systematics ( Adelberger et al. 20031). 

Returning to Fig. [21] the high-quality HIRES spectra have 
allowed us to extend the Wo vs. b relation for Lya to see 
the point where the line strength decreases rapidly, evidently 
b ~ 250 kpc. It seems likely that the rapid fall-off in Wo indi- 
cates that Lya is becoming optically thin. Indeed, under this 
assumption, the line strength in the bin with the largest b is 
Wb = 0.31 ±0.05 A which would correspond roughly to to 
log N(H I) ~ 13.25. Thus, as for the metal line transitions, a 

Lya forest decrement at z ~ 2.2. 



combination of decreasing covering fraction and linear curve 
of growth effects steepens the Wo vs b relation. 

We now explore the possibility that a simple schematic out- 
flow model might reasonably reproduce the observations of 
both the b = constraints (line strength and line shape in 
galaxy spectra) and the observed variation of absorption line 
strength with impact parameter. 

7. A SIMPLE MODEL FOR OUTFLOWS AND CIRCUM-GALACTIC 

GAS 

7.1. Line Strength and the Spatial Extent of Circum- galactic 

Gas 

Assuming a spherically symmetric gas flow with (galacto- 
centric) radial velocity outward, v our (r) [note that this is not a 
"shell"], the strength of absorption lines produced by material 
along an observers line of sight at impact parameter b will de- 
pend on the line of sight component of v out and on the radial 
and velocity dependence of the covering fraction, f c (r, v). For 
the moment, we characterize the rapid fall-off at some galac- 
tocentric radius R e ff (which may differ for each ionic species) 
as an "edge" beyond which no absorption would be detected. 

Modeling the effective value of the covering fraction (and 
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FIG. 22. — Average spectrum near Lyct in regions of HIRES spectra cen- 
tered on the systemic redshifts of galaxies with QSO line of sight impact 
parameters of 120 < b < 200 kpc (Blue; (b) = 170 kpc) and 200 < b < 280 
kpc (Red; (b) = 240 kpc). Each spectrum represents an average of 21 galaxy- 
QSO sightline angular pairs. The centroid velocities are very close to to v = 0, 
indicating that our estimates of the galaxy systemic redshifts are accurate. 
The equivalent widths fall on an extrapolation to larger b of the relationship 
obtained from galaxy pairs. 
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FIG. 23. — Diagram illustrating the simple model for estimating the relative 
interstellar absorption line strength Wo as a function of impact parameter b 
used to produce the model curves in Fig. |2T] R e ff is the characteristic size of 
the gas distribution producing saturated lines of a particular transition. Any 
systematic dependence of v ollt on radius has not been included in this model, 
but the covering fraction f c of gas giving rise to the transition of interest 
is assumed to have a radial dependence f c (r) oc r~ 7 (see text for additional 
details). 

the resulting line equivalent width) for gas along an observed 
line of sight must include an assumption about the radial de- 
pendence of the outflow velocity of absorbing material VoutM 
and f c (r, v), where the radial coordinate r takes on values be- 
tween b and R e ff, as illustrated in Fig [23] For a saturated tran- 



TABLE 5 
Wq vs. b Model Parameters 3 



Line 


7 b 


Reff (kpc) 


V ut 


fc,max 


Ly«(1216) 


0.37 


250 


820 


0.80 


C IV(1549) 


0.23 


80 


800 


0.35/0.25 d 


C 11(1334) 


0.35 


90 


650 


0.52 


Si 11(1526) 


0.60 


70 


750 


0.40 


Si IV(1393) 


0.60 


80 


820 


0.33 



a Para mete rs used to produce the model curves shown 
in Fig.|2l] 

b Power law exponent in the expression f c {r) = 

fc : niaxf 7 

c Maximum value of the covering fraction for each 
transition, measured from the composite spectrum (see 
Fig.[8) 

d Includes contributions from C IV A 1548 and C IV 
A1550 of 0.35 and 0.25, respectively. 

sition, the observed line profile (normalized intensity / vs. ve- 
locity) is given by 

/(v)=l- J f c (r,v)dl , (8) 

where v is the line-of-sight component of v out (r) relative to the 
galaxy systemic redshift, and the integral is evaluated along 
the line of sight (shaded region in Fig. [23]) If the form of 
Vout(r) is assumed, then one can associate a line of sight com- 
ponent of velocity (v) with each value of r along the sightline, 
so that f c (r, v) becomes a function of only one variable (either 
r or v). The line strength is then given by the integral of /(v) 
evaluated between the minimum and maximum line-of-sight 
velocity expected. For a velocity field which has v oul increas- 
ing with r, the extreme velocities expected would be 

I Av max(h) \ = v out (R,«)[l-b 2 /Rl s ] 1/2 (9) 

and the associated equivalent width (with continuum normal- 
ized to unity) is just 

W (v) = 2 f c (v)dv (10) 

Jo 

which can be integrated over wavelength to give the expected 
absorption line equivalent width. Assuming for the moment 
that f c (r) = fcmaxf' 1 , and that v out (r) is constant, one can ad- 
just the values of v out , 7, and 7? e ff in order to reproduce the 
observed Wo(b) for each species. The maximum covering 
fraction f CMax can be measured from the stacked spectra of 
the foreground galaxies- it is just the maximum depth of the 
absorption lines relative to the continuum for each species. 
These have been measured from an average of all 3 fore- 
ground galaxy samples, Pl(fg), P2(fg), and P3(fg), and col- 
lected in Table [5] Table [5] also contains the parameter values 
used to produce the model curves shown in Fig. [21] for each 
species. 

The examples plotted using dashed lines in Fig [21] show 
that the model does a reasonable job reproducing the rela- 
tively shallow dependence of Wo on b , and the steep decline 
at b ~ 70-250 kpc (depending on the transition). The best 
values of 7 are in the range 0.2 < 7 < 0.6, depending on the 
transition (Table [5]); note that the assumption of f c (r) oc r" 2 , 
which would apply if the absorbing clouds retained the same 
characteristic physical size as they move to larger r, appears 
strongly ruled out by the relatively shallow observed depen- 
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dence of Wo on b . Obviously, a distinct "edge" to the gas 
distribution at some radius R e ff is not physical; however, this 
radius could (e.g.) mark the point at which a transition is be- 
coming optically thin, so that both declining covering fraction 
and decreasing optical depth lead to a rapid diminution of Wo 
with b. This explanation could also account for the fact that 
Lya does not exhibit a break in the near-power-law depen- 
dence of Wb on b until b ~ 250 kpc: Lya remains strongly 
saturated as long as log N(H I) ^ 14.5, so that geometric dilu- 
tion dominates over changes in optical depth to radii beyond 
the limit of sensitivity for the the galaxy-galaxy pair samples; 
even with unity covering fraction a rest equivalent width of 
~ 0.9 A would correspond to a saturated Lya line. 

For this particular model in which v ou , is independent of 
galactocentric radius r, the combination of the measured 
fanax and v oul control the overall normalization, while the 
slope of the Wo vs. b function is dictated entirely by the 
choice of 7. As shown in Table [5] Lya, C IV, and Si IV are 
best matched if v ou , = 800 km s" 1 , while Si II and C II suggest 
v ur = 750 km s" 1 and v ou , = 650 km s" 1 , respectively. The fact 
that the latter two values are somewhat lower may be related 
to the influence of (non-outflowing) absorption near v = on 
the measurement of f c . max (see discussion in §HJ; for example, 
changing the value of f Cl max for C II from 0.52 to 0.42 (i.e., by 
less than 20%) favors v out = 800 km s" 1 rather than v out = 650 
km s" 1 . 

So far, the model has made a number of idealized assump- 
tions: spherical symmetry, uniform outflowing v eloci ty field 
(we will assess the validity of this assumption in § !7.31 l. radial 
dependence of the covering fraction f c (r) oc r" 7 , and saturated 
lines whose equivalent widths are controlled by the covering 
fraction coupled with the range of velocities expected to be 
sampled along the line of sight. It is worth a cautionary re- 
minder at this point that Wo{b) for a particular transition can- 
not be converted unambiguously into a map of total column 
density vs. b so long as the lines remain optically thick. In- 
stead, one measures the average covering fraction of gas host- 
ing that particular ion; conversely, when the optical depths 
become of order unity it becomes feasible to measure total 
column density, but at that point the constraints on covering 
fraction are diminished. 

As mentioned above, C IV, with 7 ~ 0.20, has a much more 
gradual decline with b than Si II and Si IV. This behavior 
might be attributed to a changing ionization state of more dif- 
fuse gas entrained in the outward flow, so that "new" pock- 
ets of gas become more likely to absorb in Civ as the flow 
becomes more diffuse at larger radii (i.e., it is not just the 
same "clouds" with larger cross-section). The strength of C IV 
absorption often appears "de-coupled" from the strength and 
depth of the lower ionization species in galaxy spectra (see 
Shapl evet alj 12003). which may be related to the very dif- 
ferent radial distribution of absorbing gas inferred above, as 
well as to the much larger physical coherence scale of C IV- 
absorbing gas compared to gas giving rise to lower ioniza- 
tion species observed in high resolution studies of very small- 
separation lines of sight. Nevertheless, Civ absorption usu- 
ally spans a very similar range of blue-shifted velocities in 
the galaxy spectra, strongly suggesting it is being carried in 

20 It turns out (perhaps counter-intuitivel y) th at the power law slopes re- 
lating Wb and b in the model curves in Fig. 12 II are very close to the power 
law index 7 assumed in the expression for f c (r) because of the relationship 
between impact parameter b and the line-of-sight dependence of the sampled 
velocity range. 



the same outflow. The quality of the present data does not 
justify a more detailed analysis. 

7.2. Total Absorption Cross-Section 

A very simple (and model-independent) calculation that can 
be made using the information collected in Fig.|2T|is the ex- 
pected incidence of absorption lines exceeding a particular 
threshold in Wo associated with extended gas around galax- 
ies similar to those in the current sample. This exercise is 
the converse of one undertaken many times in the past for 
various types of QSO absorbers; typically, one takes an ob- 
served value of dN/dz, the number of systems observed per 
unit redshift, and calculates the galaxy cross-section and num- 
ber density required to account for the rate of incidence, as- 
suming the cross-section is all contributed by galaxies. In 
the present case, we have a well-defined population of galax- 
ies, uniformly se lected, with a wel l -estab lished far-UV lumi- 
nosity function. Reddy & Steidel (2009) have produced the 
most up-to-date LF, determined from the same data set and 
using the same selection criteria as the galaxies used in the 
galaxy-galaxy pair analysis. They found that a Schechter 
function with 0* = 2.75 ±0.54 x 10" 3 Mpc" 3 , M* = -20.7, 
and a = -1.73 ± 0.07 is a good description of the rest-frame 
1700 A luminosity function of UV-selected galaxies in the 
redshift range 1.9 < z < 2.7. Integrating this function only 
over the luminosity range of the actual spectroscopic sam- 
ple (i.e., to 1Z < 25.5, which corresponds to L ~ 0.3L* at 
(z) = 2.3,) the number density is n ga i ~ 3.7 x 10~ 3 Mpc" 3 . Tak- 
ing b ~ 90 kpc as the point where W (C IV A 1548) ~ 0.15 A 
(see Fig. I2TI) . at the mean redshift of the foreground galax- 
ies, the co-moving cross-section for absorption per galaxy is 
(jgai ~ 0.26 Mpc 2 . The relevant co-moving path length per 
unit redshift at (z) = 2.20 is / ~ 1320 Mpc, so that the ex- 
pected number of W)(C IV)> 0. 15 A absorption systems 21 per 
unit redshift is 

^(CIV, pred) ~ (3.7 x 10 -3 ) /cr Ra i ~ 1.1 . (11) 
dz 

The observed incidence of C IV systems at the same equiv- 
alent width threshold at z ~ 2.2 is dN/dz(obs) = 2.44 ±0.29 
dSargent et alj|1988t lSteidellll990al) . so that galaxies in the 
current spectroscopic sample alone account for ^45% of 
strong Civ systems at z ~ 2. Although the statistics for 
the incidence of other lines are not as well-established, 
ISteidel & Sargent (fT992h found that, for the C II A1334 ab- 
sorption with Wo > 0.15 A and (z) = 2.35, 

^(CII, obs)~ 0.94 ±0.33. (12) 
dz 

At face value, since the extent of C II is approximately the 
same as C IV (Table [5] and Fig. |2H . comparison with eq. Q~2] 
suggests that essentially all C II A 1334 absorption with Wo > 
0.15 A is within ~ 90 kpc (physical) of a galaxy similar to 
those in the sample used here. The presence of strong C II ab- 
sorption generally depends on the self-shie lding provided by 
H I with log N > 17.2 (e.g.. lSteidellll990bl) . so that although 
we cannot measure N(H I) from the galaxy pair spectra, one 
might expect a similar cross-section for Lyman Limit System 

21 This equivalent width threshold is roughly equivalent to systems having 
log N(CIV) > 13.5 cm" 2 
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(LLS) absorption. According to lSteidelldl992lh 

-^(LLS, obs)~ 1.4 (13) 
dz 

at z — 2, suggesting that the total absorption cross-section of 
these galaxies accounts for as much as ~ 70% of all LLS ab- 
sorption at similar redshifts. 

It is important to emphasize that the absorption line statis- 
tics obtained using galaxies as background sources should 
be identical to what would be obtained using QSOs, in the 
limit of a suitably large number of QSO sightlines. In other 
words, a given average rest-frame equivalent width would be 
obtained at the same characteristic impact parameter from a 
galaxy, independent of the morphology of the background 
probes. It follows from the above discussion that most of 
the cross-section for strong low-ionization metal-line absorp- 
tion (and by extension, LLSs) is provided by rapidly star- 
forming galaxies, at least at z — 2.2. A very similar connec- 
tion between strong QSO absorbers and rapid star formation 
has been observed recently by Me nard et alj (f2009). who used 
a very large sample of Wo > 0.7A Mgll absorbers (in SDSS 
QSO spectra) over the redshift range 0.4 z 1.2, and com- 
pared with the stacked [Oil] A3727 emission signature de- 
tected within the SDSS fibers. These authors conclude that 
the strong absorption is tightly correlated with the presence 
of a rapidly star forming galaxy within ~ 50 kpc of the QSO 
sightline. By way of comparison, the statistical impact param- 
eter for producing C II absorption with Wo > 0.7 A at (z) =2.2 
is b < 60 kpc. 

We now turn to additional consistency checks by combin- 
ing the spatial information from the preceding discussion with 
information extracted from the profiles of the IS absorption 
lines in the b = spectra of typical LBGs. 

7.3. Constraints from IS Absorption Line Profiles 

From the previous section, the galaxy pair data indicate a 
covering fraction of absorbing gas with a radial dependence 
f c (r) oc r~ 7 , with 0.2 7 0.6. All measured lines except 
Lya have dropped below the current detection limit for Wo by 
b ~ 100 kpc; for larger impact parameters, spectra of higher 
resolution and S/N are required- these measurements are be- 
ing made using the bright QSO sightlines in each of the survey 
fields (e.g., Rakic et al 2010; Rudie et al 2010, in prep.). If 
we focus for the moment on scales r < 125 kpc, the dominant 
importance of the covering fraction in modulating absorption 
line strength means that there are only weak constraints on 
changes in total column density with radius. However, the 
kinematics of the absorption lines observed along sightlines 
at b = 0- the spectra of the galaxies themselves- provide a 
means of checking the models for consistency. Ideally, one 
should be able to reproduce the absorption line shape (i.e., 
covering fraction as a function of velocity) using a model for 
VoutM and f c (r) consistent with the observations of both direct 
and offset lines of sight. 

The boundary conditions provided by the direct line of 
sight (b = 0) observations are that the maximum observed 
I Avis I should be ~ 800 km s" 1 , the maximum covering frac- 
tion is typically smaller than unity, and the line shapes gen- 
erally imply that larger blue-shifted velocities are associ- 
ated with gas having smaller covering fraction. As we have 
shown, the simplest (although probably not unique) physical 
explanation would be that the higher velocity gas is located 



predominantly at larger radii 22 (see also IWeiner et al.l l2009t 
Martin & Bouche 2009). Since the outflows are believed to 
be driven by energy and/or momentum originating in the in- 
ner regions of the galaxy, a plausible model might have the 
acceleration experienced by a gas cloud depend on galacto- 
centric distance r. The details of this dependence are likely 
to be complex, since the force experienced by cool gas clouds 
would depend on the energy and/or momentum acquired by 
the gas, the ambient pressure gradient, the radial dependence 
of the clouds' physical sizes 23 , and the entrainment of ambi- 
ent gas in the flow (sometimes referred to as "mass-loading") 
as it moves outward. 

In an idealized picture of a multi-phase outflow, gas clouds, 
perhaps pressure confined by a hotter, more diffuse gaseous 
medium, are accelerated by thermal and/or radiation pressure 
to larger radii and higher velocities. The fact that the inferred 
radial dependence of the absorbing gas covering fraction re- 
mains at f c (r) oc r" 7 with 7 < 0.6 (at least for radii r < R e fi) 
implies that the clouds must expand as they move out (oth- 
erwise / c (r) oc r" 2 ). Even the most steeply declining species, 
with 7 ~ 0.6, requires that the total cross-section for absorp- 
tion per cloud increases as a c oc r 4 / 3 , or that the character- 
istic cloud size R c (r) oc r 2 ' 3 . For clouds of constant mass 
M c = 4/37rRlp c and assumed pressure gradient p(r) oc r" 2 , 
maintaining local pressure equilibrium leads to the same de- 
pendence, R c oc r 2//3 ; thus the behavior of f c (r) for Sill and 
SilV is very close to what one might expect for cooler, denser 
"cloudlets" confined by a hotter medium whose pressure de- 
creases with galactocentric radius. One might think of the 
clouds giving rise to low-ionization absorption as "expanding 
bullets", surviving as long as there is sufficient ambient pres- 
sure to confine them, or until they expand to the point that 
they no longer produce significant low-ionization absorption 
(see lSchave et alj|2007l for additional observational evidence 
of such phenomena). 

Strickland (2009) has suggested that the force experienced 
by cool clouds in a wind model obeying local pressure equi- 
librium would have radial dependence F(r) oc r" 4 / 3 ; this is 
precisely what would be expected under the circumstances 
described above (i.e., R c oc r 2//3 and p(r) oc r" 2 ). For clouds 
of constant mass, the radial dependence of cloud acceleration 
a(r) would have the same form as F{r); if additional mass 
(e.g. ambient ISM) is entrained in the wind, then a(r) would 
decrease more rapidly with increasing radius. Once the ra- 
dial dependence of cloud acceleration is specified, one obtains 
v ui(r), and thus / c (v), the expected velocity dependence of the 
covering fraction. 

Rather than develop a more detailed model for the physics 
of the multiphase outflowing material (which is beyond the 
scope of the current work), for simplicity we parametrize the 
cloud acceleration a(r) as 

a(r)=Ar~ a (14) 

where A is a constant, and seek the value of a that can repro- 
duce the IS line profiles given f c (r) inferred from the galaxy- 
galaxy pair data (Fig.|2T]and Table|5]l. Using this parametriza- 

22 Note that the simple models presented in 0above assumed a monotonic 
relationship between physical location and velocity to simultaneously explain 
the kinematics of observed Ly« emission and IS absorption 

23 Recall that the radial dependence of the covering fraction is inconsis- 
tent with clouds maintaining the same size; we show below that the cloud 
dimensions R c must increase with r. 
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tion of a(r), 

a(r) = Ar = — - — 
at 



dv{r) 



dr 



|)-*)(5l-aj) 



Assuming the clouds are "launched" from r = r m j n , equa- 
tion Q3] can be rearranged and integrated to provide an ex- 
pression for v out (r): 



VouM ■■ 



2A 
a-1 



( r l-a_ r l-a) - 5 
V mm ' ) 



(16) 



The constant A can be obtained by using the boundary con- 
dition provided by |v ma , : |, which is the velocity at which f c (r) 
approaches zero 24 and the value of R e jf estimated for the tran- 
sition of interest (Table |5j. 

We are now in a position to attempt to match observed IS 
absorption line profiles while simultaneously adhering to the 
observationally-inferred form of f c {r). For simplicity, we as- 
sume that the spectra of galaxies sample their own outflow 
kinematics at b = (i.e., we ignore the spatial extent of the 
galaxy continuum and treat observed spectra as a spatial aver- 
age over a region typically 2-3 kpc in diameter). 

For demonstration purposes, to the outflow profile is added 
an absorption component at v = with a maximum covering 
fraction varied to reproduce the v > portion of the line pro- 
file. The width assigned to this component is taken to be the 
same as the nebular line width associated with the H II regions 
in a typical galaxy (er ~ 100 km s -1 .) Fig. l24l shows a com- 
parison of the predicted velocity profiles for an IS line given 
the above model (assuming that r m j n = 1 kpc), compared to the 
observed velocity profile of the C II 1334.53 line in the com- 
posite spectrum of z ~ 2-3 LBGs with R ~ 1330. Also shown 
in Fig. [24] are the resulting relationships between r and v and 
between f c and v. The values of 7 were chosen based on the 
models (Table [5) that adequately describe how Wo varies with 
b for the relevant ion: a self-consistent kinematic model may 
be obtained for 7 = 0.35 (the value inferred for C II A1334; 
see Table [5) if a = 1.15. To provide some intuition on how the 
line profile would be altered with different parameter values, 
the predicted profile with a held fixed and 7 = 0.20 instead of 
7 = 0.35 is also plotted. It shows that assuming an incorrect 
value for the radial dependence of f c leads (in this case) to an 
absorption feature that has too much apparent optical depth at 
high values of |v| and therefore significantly over-predicts the 
values of | (v oul ) | and W). 

By altering the maximum covering fraction of the outflow- 
ing gas and the optical depth and covering fraction of the v = 
component, it is possible to reproduce a wide range of ob- 
served line profiles. Fig. [25] (dark blue model profile) shows 
another example with a chosen to produce a very strong line 
with a large blue-shifted centroid such as those observed in 
the sp ectrum of the lensed LBG MS1512-cB58 dPettini et al.l 
2002). In this case, it was necessary to change f c ,max to unity 
and to adjust the function a(r) to be steeper than for the ex- 
ample in Fig. [24] Since the line is again C II A1334.53, once 
again 7 = 0.35. The orange model profile is that obtained 
when all parameters are held fixed but a is reduced to the 
value that provided a good fit to the profile in Fig. [24] Qualita- 
tively, as the function a{r) becomes steeper, the covering frac- 
tion at high velocity increases along with the line equivalent 

24 In other words, this is the maximum velocity beyond which the covering 
fraction is too small to produce absorption that can be distinguished from the 



width (given the imposed boundary conditions); the physical 
meaning of a steeper a(r) is that material is accelerated more 
rapidly, such that it reaches high velocity before it becomes 
geometrically diluted. 

The bottom panels of Figs. l24l and |251 are also interesting to 
consider. They show that, at least in the context of the cur- 
rent models, the outflowing gas attains close to its maximum 
velocity within ~ 3 — 10 kpc (depending on the model) of the 
center of the galaxy, where whatever mechanism is at work in 
accelerating the clouds is most effective. At larger radii, the 
outflow velocities change slowly with increasing radius. Re- 
call that our simple model for inferring f c (r) from Wo vs. b 
assumed that v out is independent of r. Since the current model 
requires a(r) = dv ou ,(r)/dt oc r~ a with 1.15 £ a ^ 1.95, most 
of the acceleration up to the maximum observed velocity oc- 
curs when r « 10 kpc . Th us, we conclude that having as- 
sumed a fixed v ou , as in jj7.1l was justified, given other sources 
of uncertainty. 

We caution that in the particular family of models described 
above, the parameters 7 and a are covariant and combinations 
of the two can usually be found that will adequately reproduce 
observed IS line profiles. It is only because we have informa- 
tion on the radial dependence of f(r) from the galaxy-galaxy 
pairs (constraining 7) that anything like a unique solution re- 
sults. Since the value of 7 is statistical (and so may not apply 
for particular individual galaxies) it would be unwise to read 
too much physical significance into the inferred values of a 
needed to fit a particular line profile. For example, the C II 
line profile in Fig.[24]can be well-fit a model with the param- 
eters (a, 7) = (1.33,0.60), and that shown in Fig. [25] can be 
equally well-fit by the combination (0,7) = (1.25,0.20). 

8. DISCUSSION 
8.1. Connection Between b = and b » Sightlines 

The models discussed above are schematic only, and prob- 
ably not unique in adequately describing the relatively crude 
observations. Nevertheless, we have seen that a reasonably 
successful model for the kinematics and overall geometry of 
the outflows includes a characteristic asymptotic velocity of 
~ 700 - 800 km s" 1 . Such high velocities are required to ex- 
plain the shapes of IS absorption line profiles "down the bar- 
rel" for b = sightlines, as well as the strength of the IS ab- 
sorption and its relatively shallow dependence on impact pa- 
rameter for sightlines with b > 10 kpc. This suggests that the 
absorption at large b is causally related to the outflowing ma- 
terial observed at b = 0, in which case the gas at large r would 
most naturally be identified with cool gas carried out in an ear- 
lier phase of the current episode of star formation. If this were 
the case, one might expect the geometry of the metal-enriched 
gaseous envelope of a galaxy to be time-dependent. 

For example, galaxies with spectral morphology similar to 
that of MS1512-cB58 (and represented by the model shown 
in Fig. [25l l tend to be those with estimated stel lar population 
ages (star formation episodes) of t s f < lOOMyr ( Shapl ey et alJ 
l200UlKornei et alj2009l) . Even for gas moving at 800 km s _1 , 
reaching r ~ 100 kpc requires ~ 120 Myr, i.e. longer than 
the inferred duration of the current star formation episode 
and ~ 20% ofUf typ ical of UV-selected galaxies at z — 2-3 
(Shapl ey et alJ 12005). Since the inferences about the CGM 
geometry from the galaxy-galaxy pairs are statistical, and 
the subset of the pair sample having ancillary information 
on stellar populations is small, we are not in a position at 
present for a definitive test. Significant variations in star for- 
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FIG. 24. — (Top): Model line profiles (convolved with the appropriate instrumental resolution) for an IS absorption line with radial covering fraction f c (r) oc r 7 , 
maximum covering fraction for outflowing material fc.max = 0.59, and radial cloud acceleration a(r) oc r~ a , compared to the C II A1334.53 profile from a composite 
of R = 1330 galaxy spectra. Two models are plotted, both of which include an absorption component with v = 0, <r v = 80 km s _1 , and fc.max = 0.4 : the dark 
blue profile assumes 7 = 0.35 as inferred from the galaxy pair measurements for the C II line, and a = 1.15, which provides a good match to the observed line 
shape. The orange curve represents the model line profile if 7 = 0.2 and all other parameters are held fixed. Note that the relatively small change in 7 produces 
a line profile inconsistent with the data, with a high velocity wing that is too strong. Both are constrained to have v max = —700 km s -1 at the point where the 
line profile becomes indistinguishable from the continuum. The covering fraction of outflowing material is assumed to go to zero at r = R e /f = 90 kpc Bottom: 
The dependence of the covering fraction (purple, righthand axis) f c (v) and the relation between v m and galactocentric radius r (cyan, left-hand axis). The solid 
curves correspond to the dark blue profile, with the dashed curves indicating values for the orange profile. Note that the velocity profiles are identical for the 
two models, since they use the same value of a. The red curve shows the estimated escape velocity as a function of galactocentric radius r for an NFW halo of 
mass 9 X 10" Mq ; the green dotted line is the virial radius r,„, for the same halo. Note that at r = r,„> (in the context of this model) the outflow velocity exceeds 
Vescvvir). 



mation rate or of the mass flux of outflowing material might 
also be expected; if star formation were to shut down com- 
pletely, the corresponding diminution of outflowing gas might 
eventually be recognized as "gaps" in the CGM near to the 
galaxy, but lingering high-velocity material at larger radii. 
Such remnant wind material may have been observed in the 
spectr a of massive post-sta rburst galaxies at intermediate red- 
shift dTremonti et al.l 120071) . At z — 2-3, our sample does 
not currently include galaxies without significant star forma- 
tion, though it does include some galaxies with small inferred 
gas fractions and low specific star formation rate SFR/M* 
dErb et ai1l2006bl) . Unfortunately, these too do not have high 
enough surface density to have contributed significantly to our 
galaxy-galaxy pair samples. 

One of the expectations of the simple kinematic model out- 
lined above is that most of the absorption seen in the spectrum 
of a LBG is due to gas located within ~ 10 kpc, while most 
of the absorption observed in offset lines of sight is due to gas 



that makes a relatively small contribution to the line profile in 
the b = spectra. At large r, because gas giving rise to a par- 
ticular absorption line transition has lower f c , the constraints 
on its velocity from the b = line profiles become progres- 
sively weaker. This means that there could be lower or higher 
velocity material at large radii and small covering fraction that 
might go unnoticed compared to material with the same veloc- 
ity but much smaller galactocentric radius (because it would 
have a less-diluted f c .) 

In this context, it is interesting to consider what is responsi- 
ble for the "extra" IS absorption in the higher baryonic mass 
sub-sample discussed in §|4] above. It is intriguing that the 
strength of the v ~ components of IS absorption is such 
a strong function of galaxy m ass in both the curr ent sample 
and in the z ~ 1.4 sample of Weineret al. (2009p. A com- 
pari son of the highest an d lowest stellar mass bins in Table 
1 of IWeiner et al.l (12009) shows that the line strength of the 
symmetric component increases with stellar mass by a factor 
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FIG. 25. — Same as Fig. [24] but with increased maximum covering fraction of the outflowing gas f c ,max = 1, and a slightly larger covering fraction f c ,max = 0.5 
for the component at v = 0. The dark bl ue m odel profi le shows th at for 7 = 0.35, a good match to the observed C II line in the spectrum of the lensed LBG 
MS1512-cB58 (R ~ 1500, to match Fig. 1241 see Pettini et al. 2000) is obtained if a = 1.95. The orange curve illustrates the change in the model line profile if 
a = 1.15 is assumed i nstea d, with all other parameters held fixed. Clearly the latter model fails to produce the depth of the absorption at high values of |v|. As 
for the example in Fig. 1241 the outflow velocity is significantly higher than v esc at the virial radius of the assumed halo. 



of at least 16.6. iWeiner et al.1 J2009h argue that the absorp- 
tion near v = is due primarily to stellar photospheric Mg II 
absorption. However, the set of far-UV transitions used for 
our z ~ 2.3 sample are not contaminated by significant stel- 
lar photospheric absorption, and yet the line profiles exhibit a 
very similar trend with galaxy mass. We therefore conclude 
that the "excess" absorption near v = arises from IS gas, at 
least in our z ~ 2 - 3 sample.. 

There is generally little evidence for infalling (redshifted) 
gas in the far-UV spectra of the galaxies, and even less ev- 
idence that v(r) could be decreasing with increasing r (the 
line shapes would be very different if this were the case). It 
would also be difficult to explain the spectral morphology of 
Lya emission (see § HJ if the relationship between v out and r 
were inverted, as might be expected if stalling winds and/or 
infall dominated the gas flows. The expected signature of in- 
fall as seen in the galaxy spectra is IS absorption velocities 
redshifted with respect to the galaxy systemic velocity. Such 
absorption would be expected if gas were prevented from es- 
caping the galaxy potential and eventually began to fall back 
onto the central regions; a similar signature could also result 
from "fresh" (rather than recycled) infalling cool gas that is 
expected to have p hysical conditions (T ~ 10 4 - 10 5 K; e.g., 
iGoerdt et al.1 12009) very similar to the cool outflowing gas. 



The absorption line strength for the metallic species we have 
been able to measure would be relatively insensitive to the 
metal abundance of infalling material, once again due to the 
high degree of saturation in the strong IS lines. In other 
words, the presence of infalling gas should be easily detected 
in the galaxy spectra (both at b = and at large galactocen- 
tric radii). Redshifted absorption should be observed in the 
velocity range £j v + 300 km s" 1 as long as the covering 
fraction is within a factor of a few of that inferred for outflow- 
ing gas. Under the current model, infalling gas should have 
increasingly redshifted velocities as its covering fraction in- 
creases (i.e., as r becomes smaller); the expected line profiles 
would then be inverted with respect to the observations, in the 
sense that line profiles would be deepest at the high velocity 
end of the distribution. 

We have discussed evidence from the composite spectra in 
$4]ihat galaxies with higher baryonic mass are more likely to 
have stronger absorption at v ^ 0, which could possibly be a 
signature of stalled outflowing materia l and /or infalling gas. 
This possibility is discussed further in {38.4l below. 

8.2. The CGM and Dark Matter Halos 

Fig- EH shows that among the strong lines we can measure 
from the composite galaxy spectra, all except Lya appear to 
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become much weaker at impact parameters of b >, 70- 100 
physical kpc even though the rate of decline of Wo(fo) (at 
smaller values of b) varies among the observed lines. It is 
interesting to ask whether this particular scale has physical 
significance (as opposed to being an artifact of the limited 
sensitivity of the data) given that we know something about 
the properties of the galaxies associated with the extended 
gas. The particular objects in the present sample, drawn from 
a spectroscopic survey of "BX" UV-color-selected galaxies, 
have b een well-characterized in terms of their spatial clus- 
tering (Adelberger et al. 2005b) a nd stellar population pa- 
rameter s (lErb et al.ll2006ct IShaolev et"aTl 12005b iReddy et all 
120051) . IConrov et al.l d2008l) hav e used the observational re- 
sults of Adelber ger et al.l d2005bl) to match these z — 1.9 — 2.6 
galaxies to dark mat ter halos in the Millennium simulation 
(Springel et al. 2005). They find a good match to both the 
clustering strength and the space density of the "BX" galax- 
ies with dark matter halos having Mh a i > 4.2 x 10 11 M Q and 
(Mhaio } = 9 .0 x 1 1 1 M Q . If we assume that the dark matter ha- 
los "formed" several hundred million years prior to the epoch 
of their observation (e.g., a galaxy observed at z ~ 2.4 might 
have formed at z ~ 2.8 for a typical inferred stellar popula- 
tion age of ~ 500 Myr), one can estimate the corresponding 
virial radius 25 , which is ~ 64 kpc for the minimum halo mass 
of 4.2 x 10 11 M Q and ~ 82 kpc for a halo with the expected 
average mass of 9 x 10 11 . These values of r V! > would increase 
by ~ 10-15% if the "formation" redshift were assumed to be 
Z = 2.3 instead. It is possible that the similarity of r,,„- and the 
radius at which the Wo vs. b curve begins to steepen for some 
of the low ionization absorption lines is not a coincidence. 
If the clouds are confined by the pressure of a hotter, more 
diffuse medium, then r V! > may be approximately the radius at 
which the clouds become too diffuse and highly-ionized to 
produce significant columns of low-ionization metals. 

By assuming a radial density profile for the characteristic 
dark matter halo, we can estimate the "escape velocity" v esc ('") 
from a point at radius r re lative to the center of mass. Assum- 
ing a lNavarro et al.l ([1997) (NFW) dark matter density profile, 



PDM(r) oc 



1 



'l+cr/r vir ) cr/r vir 



(17) 



the corresponding escape velocity from radius r is given by 



2GM vir ln(l + cr/r vir ) 
r ln(l+c)-c/(l+c) 



(18) 



where c is the NFW concentration parameter. Assuming a 
"formation" redshift ofz = 2.8, M vir = 4.2 x 10 11 M©, andr= 1 
kpc, Vesc — 561 km s" 1 , while for the expected average halo 
mass from above, v esc ~ 730 km s -1 , both assuming c = 7. 26 
Perhaps more relevant, the escape velocity for material at r = 
r VI > is significantly lower, v esc ~ 310 km s -1 for the minimum 
mass halo (r V! > = 64 kpc) and v esc ~ 400 km s -1 for the average 
mass halo (r,,,y = 82 kpc). The lower panels of Figs. l24l and 
[25] show the estimated v esc as a function of r; note that in the 
context of the model, both cases have v out (r vlr ) > v f , JC (r v( >). 

25 Here we are assuming that the virial radius is the radius within which 
the average halo matter density is ~ 178 times higher than the mean density 
of the universe at the time of initial collapse. 

26 Vesc is relatively insensitive to c, changing by < 10% as c goes from 5 to 
9, but is somewhat more sensitive to the value assumed for z/ (e.g., assuming 
Zf ■ = 4 changes the velocities for the same halos above to v esc (r = 1) ~ 640 
km s~ ! and 830 km s -1 respectively.) 



8.3. Comparison to Observations at Lower Redshifts 

The observed value \v max \ ~ 800 km s -1 so common among 
the galaxies in the UV selected samples discussed above sug- 
gests that much of the outflowing gas is destined to become 
unbound from the parent galaxy. One might naturally ask why 
v max seems to be so consistent within the current z ~ 2 - 3 
spectroscopic samples. For halos forming at a given red- 
shift, the expected dependence of v esc on mass is fairly shal- 
low, v esc (r V i r ) oc M[?; ( 3 . The inferred range of baryonic mass 
Mbar (stars + cold gas) among the "BX" sample at z — 2-2.6 
is reasonably described by a log-norm al distribution with 
log M har = 10.61 ±0.34 (lErb et al.li2006d) : assuming the cos- 
mic baryon-to-dark matter ratio for all of the galaxies, the ex- 
pected range in v esc among the sample would then be ±25%, 
which is marginally compatible with the observed consistency 
of v max among the BX sample of galaxies. Similarly, the range 
of star formation rates a mong the NIRSPEC Ha sample is 
SFR ~ 30 ± 15 M Q yr" 1 (lErb et al l2006bl) meaning that if the 
con-elation of v max oc SFR 25 (found by IWeiner etail (120091) 
for galaxies at z ~ 1 -4) applied to the z ~ 2 - 3 galaxies, the 
expected variation of v max would be ±6%, which is indistin- 
guishable from "constant" given the uncertainties in the mea- 
surement 

IStrickland & Heckmanl d2009l) have emphasized that out- 
flow velocities observed in cool gas entrained in super-wind- 
driven outflows will be significantly smaller than that of the 
hot "wind fluid", which for M82 they infer to have charac- 
teristic velocity of 1400-2200 km s" 1 based on hard-X-ray 
observations. The highest velocity observed in cool material 
(via Ha filaments) in M82 is v out ~ 500-600 km s _1 . The 
typical galaxy in our sample has a bolometric luminosity of 
~3x 10 11 Lp), roughly 5 times highe r than that of M82 (cf. 
IReddy et al.ll2006t iSanders et al.ll2003T) . Applying the scalin g 
expected for momentum-driven winds 27 (Mu rray et al.l 2005). 
Vmax oc L 025 , by analogy with M82 one would expect the cool 
component of the outflows for a typical z ~ 2 - 3 galaxy in 
our sample to reach 750-900 km s" 1 . This range is clearly 
consistent with the observations. 

8.4. Evidence for "Cold Accretion" ? 

It is clearly of interest to understand the expected obser- 
vational signature of the accretion of gas in the context of 
the proposed schematic model of the CGM. The current mod- 
els of cold accretion predict characteristic gas temperatures 
of - 1 -5 x 10 4 K (e.g. iGoerdt et alj [20091) . which is the 
same range of temperatures expected for the cool outflowing 
material seen in absorption against the galaxy far-UV con- 
tinuum. The bulk of the cool accreting gas is predicted to 
lie within c± 100 kpc for halos of mass 10 12 M Q . According 
to iDekel et all (2009), absorption due to cold stream accre- 
tion should have a covering fraction as seen by background 
sources of ~ 25% for 20 < r < 100 kpc with N(H I)> 10 2() 
cm" 2 , with infall velocities ~ 200 km s" 1 . Along lines of 
sight with b » 0, one might expect absorption centered 
near v/ os ~ relative to the galaxy systemic redshift (because 
on average the sightlines should intersect both blue-shifted 
and red-shifted IS material), with velocity range of perhaps 
±150 km s" 1 . Sightlines with b ~ (i.e., in the spectra of 

27 IMurrav et alj (2005) point out that cool clouds entrained in a hot flow 
behave like momentum-driven winds whether the source of the flow is pri- 
marily momentum-driven or energy-driven. 
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the galaxies themselves) should intersect inflowing gas with 
< v, OJ < + 200 km s" 1 . 

As mentioned in 34.31 above (and illustrated in Fig. ITTb . 
the "excess" absorption in the composite spectrum of galax- 
ies with M[, ar £j 4 x 10 10 M may be consistent with both the 
expected kinematics and optical depth (or covering fraction) 
for infalling gas. It is notable that the v/ 0J > component of 
absorption is seen in low-ionization species, but not in more 
highly ionized species such as C IV that are otherwise present 
over all velocities V; os < km s" 1 . This may imply that the 
gas in the V/ 0J > kinematic component is more self-shielded 
than the IS material that contributes most of the observed line 
equivalent width, suggesting it lies primarily at small r. In any 
case, even when it is present, a putative infalling IS compo- 
nent cannot account for the observation of strong absorption 
in both Lya and metal lines whose strength requires high ve- 
locity dispersion and covering fraction. Nor can accretion ex- 
plain the observed kinematics of Lya emission discussed ex- 
tensively in ij5] Nevertheless, we know that galaxies with halo 
masses of ~ 10 12 M Q must accrete material in some form over 
the typical star formation timescale of ~ 500 Myr in any rea- 
sonable hierarchical model; infall of low-metallicity gas may 
also be necessary for producing the observed mass-metallicity 
relation at z ~ 2 (e.g. lErbll2008l) . 

Possibly more intriguing than the dominance of outflowing 
IS material is the observation (bot h in the present (z) = 2.3 
sample and in the z ~ 1 .4 sample of lWeineret al. (2009])) that 
galaxies below a threshold in baryonic mass appear to lack 
the Vios > IS absorption component. In the z ~ 2.3 sample 
the mass threshold appears to be near M\, ar ~ 4 x 10 10 M Q , 
below which essentially all IS absorption is blue-shifted. We 
have seen that the probable range in dark matter halo mass 
(assuming correspondence between baryonic mass and halo 
mass) for the lower-M/ M( - half of the Ha sample is 4 x 10 11 
M© M/ m i ^,9x 10 11 Mq- a range in mass over which cold 
accretion is believed to be near its peak at z — 2. Paradox- 
ically, the sub-sample showing possible evidence for an ac- 
creting component of cool gas is the half with M/ lfl / ^ 10 12 
M Q , i.e. close to the critical mass at which virial shocks be- 
gin to prevent cold g as from streaming to the central regions 
of the galax y (e.g.,[D ekel & Birnboim 2006; ICeverino et al.l 
2009, but cf. Keres et al. 2009). Lower baryonic mass galax- 
ies with similar star formation rates show no evidence for an 
infalling component with significant covering fraction, at least 
as compared to that of out-flowing material. The extended 
nature of the galaxies in whose spectra the IS lines are mea- 
sured provides a spatially-averaged IS absorption profile to 
detect both the coolest regions of inflowing material as well 
as more highly ionized gas (i.e., having N(H I) as low as 10 14 
cm -2 ) that would present a larger covering fraction. It would 
be very hard to miss infalling cool gas, which would appear 
as redshifted absorption for b = sightlines. Moreover, sight- 
lines at b » would produce lines too weak to be consistent 
with the Wo vs. b results from §|6]due to the quieter velocit y 
fields expected for infalling material; see lDekel et al.l ([2009). 

With the possible exception of the v > absorption in the 
more massive half of the current sample (which could have 
alternative explanations; see discussion in §0, the observa- 
tions reveal an absence of evidence supporting cold flow ac- 
cretion as currently envisioned. It is possible that the streams 
cover a much smaller fraction of 4ir steradians than the 
~ 20-25% estimated fr om simulations dGoerdt et al.l 120091: 
iDekel & Birnboiml20 08) in order to remain undetected by the 



absorption probes; however, since the absorption line probes 
are sensitive to H I column densities as low as N(H I) <, 10 15 
cm" 2 , they would be sensitive to much more than the coolest, 
most highly collimated regions of the accretion flow. In any 
case, there seems to be no way to reconcile the observed CGM 
absorption line strength and kinematics with the results of 
simulations which seem consistently to predict that accretion 
of cool gas should be dominant over outflow for galaxies with 
M tot ~ 10 12 M Q . Taken at face value, it seems that the impor- 
tance of cold accretion has been significantly over-estimated 
- or at least that its observational signature must be more sub- 
tle than suggested by the early predictions. Of equal concern 
is that the influence of outflows, affecting large regions of the 
CGM of relatively massive galaxies, seems to have been se- 
riously under-estimated; in many cases, outflowing material 
has been completely ignored or deemed negligible for galax- 
ies in the same mass range as in our sample. There is strong 
empirical evidence for high velocity outflows whose influence 
extends to galactocentric radii of at least 125 kpc around an 
average galaxy in our sample- roughly the same range of r 
over which cold accretion is supposed to be most observable. 
It would clearly be interesting to understand how high veloc- 
ity outflows would interact with cool accreting gas in the case 
that both are occurring simultaneously in the same galaxies 
and where the relevant physics of both processes have been re- 
alistically modeled. At present, the bulk of the observational 
results do not support the direction in which the theory has 
been moving- a situation which clearly needs to be resolved. 

8.5. Lya Emission Re-Visited 

In $5] we presented a simple 1-D model in an effort to 
understand the kinematics of Lya emission in the galaxy 
spectra. The model assumed that the observed blue-shifted 
IS absorption is providing information on the relevant gas- 
phase kinematics on the "far side" of the galaxy. We 
showed that the scattering of Lya photons in circumgalac- 
tic gas having large bulk velocities and steep velocity gra- 
dients would produce redshifted Lya emission similar to 
that observed. A somewhat more physically-motivated CGM 
model was presented in $7] The outflowing material is 
multi-phase, where the covering fraction of gas giving rise 
to a particular transition (rather than the optical depth) ver- 
sus velocity and galactocentric radius modulates the line 
profile shapes and the line strength as a function of im- 
pact parameter. A more sophisticated model for Lya emis- 
sion, including photon di ffusion in 3-D physical space as 
well as velocity space (e.g., ILaursen & Sommer-Larsen 2007; 
IVerhamme~ et al. 2008j), would be needed for a detailed under- 
standing of the Lya morphology. Such modeling is beyond 
the scope of the present work; however, the multi-phase gas 
involved in galaxy-scale outflows is likely to lead, qualita- 
tively, to Lya kinematics similar to the observations. 

As discussed in §|5] the high velocity and (more impor- 
tantly) the large velocity range evident from the observed IS 
absorption lines allow Lya photons initially produced by re- 
combination in the galaxy H II regions to work their way out- 
ward in both real space and velocity space via multiple scat- 
tering events. Rather than a shell or a continuous medium of 
high H I optical depth, the velocity gradient due to the accel- 
eration of the clumpy outflowing gas allows Lya photons to 
migrate to larger radii and higher velocities until they are suf- 
ficiently redshifted that photons emitted toward an observer 
can escape the galaxy and its circumgalactic gas. Since we 
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know that the clouds giving rise to low-ionization absorption 
lines typically extend to at least r ~ 70 kpc, the model we have 
proposed implies that there should be a diffuse Lya halo on 
the same physical scales. The total luminosity in the diffuse 
halo may, in some cases, represent a significant fraction of 
the galaxy's production of Lya, albeit distributed over a pro- 
jected surface area as much as >, 1000 times larger than the 
galaxy continuum light. Of course, most of this Lya emis- 
sion would not be included in the narrow slit typically used 
for galaxy spectroscopy. It is possible that the diffuse emis- 
sion has already been obs erved in very deep Lya images in 
the fo rm of Lya "Blobs" (ISteidel et all 12000 : Matsu daetail 
2004)- the most extreme examples- or as diffuse emission 
eviden t only after stackin g spatial regions surrounding normal 
LBGs (Havashi no et alj |2004. Steidel et al, in preparation). 

The main point is that the structure and kinematics of the 
CGM ultimately control the morphology of escaping Lya 
emission, and the Lya "photosphere" of a galaxy is expected 
to be roughly coincident with the distribution of gas responsi- 
ble for the observed IS absorption features. This "prediction" 
of very extended Lya halos on scales of ^50-100 kpc is 
remarkably similar to recent predictions of cooling Lya emis- 
sion associated with cold accretion onto galaxie s of very sim- 
ilar mass scale to those in t he current sample (Goerdtet al. 
2009- lDiikstra & Loebl20 09). Lya emission is expected to be 
both redshifted and relatively axisymmetric for Lya originat- 
ing in the galaxy and scattering through the CGM, whereas 
Lya emission from accreting gas should favor blue-shifts and 
will be concentrated in a few dense filaments. However, it re- 
mains unclear how observationally distinct the extended Lya 
emission arising from these two very different processes will 
be. 

We intend to address these issues in greater detail in a sep- 
arate paper. 

8.6. Outflow Mass Flux and the Gas Content of the CGM 
In principle, the type of outflow model described above can 
be used to estimate the outward mass flux M ou , as well as 
the total (cool) gas content of the CGM, Mqgm- Mass out- 
flow rates o f cool gas hav e been estimated for nearby star- 
bursts (e.g., iMartinl 120051) , for rare ind ividual examples of 
high redshift galaxies (e.g., Pettini et al. 2000), and for com- 
posite spectra of star-for ming galaxies at intermediate red- 
shifts dWeiner et al.ll2009l) . In all cases, the mass flux is esti- 
mated to be M ou , ^ M* where is the galaxy star formation 
rate. The calculations generally assume that the outflowing 
material is in a thin shell located at a galactocentric radius r 
just beyond the observed stars, traveling with outward veloc- 
ity given by the centroid of the outflowing component of the 
IS absorption lines, where the total gas column densities have 
been estimated from the strength of low-ionization metal ab- 
sorption lines (for z < 1.6) or from the observed N(H I) mea- 
sured from the Lya absorption, in the case of MS1512-cB58. 
The uncertainties in these calculations are very large even in 
the context of the simple shell model, since they depend on 
the assumed shell radius as r 2 and the total gas columns are 
poorly constrained even when N(H I) has been directly mea- 
sured, due to unknown ionization corrections. 

Unfortunately, the less-idealized model introduced above 
does not necessarily improve the situation; the IS absorption 
is produced by gas with a large range of both radius and ve- 
locity, rather than a shell at a single radius. This means that a 
galaxy spectrum contains information on the integral absorp- 



tion along its line of sight, and therefore includes material 
deposited by both current and past winds. Since the depth of 
absorption lines is significant out to R e ff ~ 100 kpc in our 
outflow model, gas that was launched up to £ 100 Myr ago 
may still contribute to the line profiles, albeit diluted by a de- 
creasing covering fraction with increasing r. 

In the context of our model, cool gas deposited by outflows 
to r ~ 100 kpc represents an appreciable fraction of the ex- 
pected Mbar.tot — 1 .5 x 10 11 M Q associated with a typical dark 
matter halo mass of ~ 9 x 10 11 M . If we assume for def- 
initeness that R e ff ~ 80 kpc and that the asymptotic outflow 
velocity is v out ~ 750 km s" 1 , a packet of cool gas launched 
from r,„,„ ~ 1 -2 kpc will reach r ~ 80 kpc after % ~ 1 -3 x 10 s 
years. If there is a steady-state flow with mass flux M out , the 
total cool gas content within R e fj = 80 kpc deposited by the 
outflow would be given approximately by 



(M CC m) = {M ou t)tw ■ 



(19) 



The model that best reproduces the radial behavior of f c for 
low-ionization species in the CGM of typical galaxies implies 
cloud densities that increase with galactocentric radius as p c oc 



r ( 137.31 ). If we assume that the initial number density of 
atoms in a typical cloud is no cm" 3 at the base of the flow r„, 



then the total cool gas mass between r„ 
kpc is 



1 kpc and r = 80 



M CG u(r < 80) ~ 3 x 10 10 n M , 
(M CC m) 
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Note that this estimate does not attempt to account for hot 
or highly ionized gas that is also likely to be associated 
with supernova-driven outflows (e.g., IStrickland & Heckmanl 
120091) . 

An estimate of «o follows from our earlier assertion that 
R e jf for low-ion species likely corresponds to a threshold 
N(H I) XL 10 17 cm" 2 (i.e., a Lyman Limit system). The neutral 
fraction for such H I column density gas ionized by the meta- 
galactic radiation field is estimated to be ~ 10~ 3 (e.g.,[Steidel 
1990b), so that the total H column will be ~ 10 20 cm" 2 . For 
a line of sight with b ~ 80 kpc, the corresponding average 
density «//(80) ~4x 10" 4 cm" 3 . In the case of the aforemen- 
tioned r" 2 density dependence, 



R. 



eff 



-, 2 



n H (80) - 1 cm" 



(22) 



Adopting «o = 1 cm" 3 , cool gas in the CGM would account 
for at least 20% of the total baryons associated with a typical 
galaxy, with mass comparable to the total Mt, ar (cold gas plus 
stars) in th e central few kpc , which has (Mt, ar ) ~ 4 x 10 I() M 
at z = 2.3 (Erb et al. 2006c). If one assumes that the typical 
stellar mass M* ~ 2 x 10 10 M Q was formed over the same 
timescale (i.e., f 80 ~ 1 .3 x 10 8 yr), then (SFR) ~ 150 M yr" 1 . 
Thus, a crude (but largely independent) calculation yields the 
same qualitative result obtained with shell models in previous 
work: i.e.,M „, ^M*. 

In any case, it is difficult to escape the conclusion that the 
CGM contains a substantial fraction of galactic baryons at 
z ~ 2 - 3, and that the gas-phase kinematics and geometry 
points to a causal connection between rapid star formation 



The Circum-Galactic Medium at z = 2 — 3 



31 



and the presence of large amounts of gas at large galactocen- 
tric radius 28 . 

9. SUMMARY AND CONCLUSIONS 

We have used a relatively large sample of 1.9 ^Sz^ 2.6 
galaxies with accurate measurements of systemic redshift z svs 
and reasonably high quality rest-frame far-UV spectra to ex- 
amine the relationship between galaxy properties and UV 
spectral morphology. Our main focus has been on the kine- 
matics and strength of interstellar absorption and Lya emis- 
sion and their implications for the galaxy-scale outflows ob- 
served in all rapidly star-forming galaxies at high redshifts. 
Using this well-observed subset as a calibration, we then com- 
bined the rest far-UV galaxy spectra drawn from a much- 
larger parent sample with additional spatial information for 
the same ensemble of galaxies provided by the spectra of 
background galaxies with small angular separations. Using 
these joint constraints, we constructed simple models of the 
kinematics and geometry of the "circumgalactic medium"- 
the interface between star forming galaxies and the IGM. Our 
principal conclusions are as follows: 

1 . We have used the Ha sample of z — 2 - 2.6 galaxies to- 
gether with their far-UV spectra to produce a revised calibra- 
tion that allows deriving systemic redshifts for star-forming 
galaxies from their far-UV spectral features (strong IS lines 
and Lya emission). In the absence of stellar photospheric ab- 
sorption lines or nebular emission lines in the rest-frame op- 
tical, the most accurate estimates of the systemic redshift are 
derived from the centroids of strong IS absorption lines. Ap- 
plying a shift of Av = +165 km s" 1 (Az = +0.0018) at z = 2.3 
to a measured absorption redshift provides an estimate of a 
galaxy's systemic redshift accurate to ±125 km s" 1 ; redshifts 
measured from the Ha emission line (when available) have a 
precision of ±60 km s" 1 . Both methods for estimating red- 
shifts have no significant systematic offset relative to the red- 
shift defined by stellar (photospheric) absorption features. 

2. Using only the Ha sample of BX galaxies, we found 
mean velocity offsets for the centroids of the strong IS lines 
and Lya emission line of Avis = -164 ± 16 km s _1 and 
Av LyQ = +445 ± 27 km s" 1 , respectively. We searched for 
significant correlations between the kinematics defined by the 
centroid velocities of IS and Lya lines and other measured or 
inferred galaxy properties. Within the Ha sample, the only 
significant correlations found were between Avis and galaxy 
mass estimated from the sum of inferred stellar and gas mass 
(Mhar) as we ll as the independently estimated Mj yn . The sense 
of the correlation is that Avis is smaller (i.e., less blue-shifted) 
in galaxies with larger masses. 

3. Despite the trend described in point 2 above, the velocity 
\vmax\ of the maximum blue-shift observed in the absorption 
profiles is essentially identical for the sub-samples of galax- 
ies with M bar > 3.7 x 10 10 M and M bar < 3.7 x 10 10 M , 
with \v max \ ~ 800 km s" 1 . However, the higher-mass subset 
has both Avis and Av LyQ shifted toward positive velocities 
(i.e., more redshifted) by ~ 200 km s" 1 . The differences in 
the composite IS and Lya profiles can be explained by an 
additional component of absorption at v ^ that appears to 
be absent in the lower-mass sub-sample. The extra absorp- 
tion, which has a peak near v = and a centroid velocity of 
v ~ +150 km s" 1 , could conceivably be a signature of infalling 
gas or stalled winds falling back on the galaxy, but is most 

28 See|Menard et al. 2009 for additional support for this causal connection. 



likely to be gas at small galactocentric radii based on the line 
profiles. In any case, the velocity centroids of the IS lines, 
which are commonly used as a proxy for the "wind velocity" 
for high redshift galaxy samples, are actually modulated al- 
most entirely by gas that is not outflowing. 

4. The Lya emission profile, like (Avis), is modulated by 
the covering fraction (or apparent optical depth) of material 
near v ~ 0. We show that, by using the information on the 
covering fraction near v ~ 0, the kinematics of blue-shifted 
gas, and assuming spherical symmetry, the behavior of Lya 
emission with respect to IS absorption can be reproduced us- 
ing simple models. In general, the red wing of Lya emission 
is produced by Lya photons scattered from outflowing gas on 
the opposite side of the galaxy. The apparent redshift of the 
Lya centroid is a manifestation of the fact that only photons 
scattering from material having a (redshifted) velocity large 
enough to take the photons off the Lya resonance for any ma- 
terial between the last scattering and the observer can escape 
in the observer's direction. The spectral morphologies of Lya 
emission and IS absorption are most easily understood if the 
highest velocity material (either redshifted or blue-shifted) is 
located at the largest distances from the galaxy, i.e., that v(r) 
is monotonically increasing with r in a (roughly) spherically 
symmetric radial flow. The same geometric picture works 
well to reproduce both the line-of-sight (b = 0) absorption pro- 
files and the absorption line strength as a function of impact 
parameter (see point 6 below). 

5. We have demonstrated that the use of the far-UV spec- 
tra of galaxies within projected angular pairs (with discrepant 
redshifts) provides an opportunity to constrain the physical 
location of circumgalactic gas around typical galaxies. A set 
of 512 such galaxy pairs, on angular scales ranging from 1" 
to 15" (~ 8 - 125 kpc), have been combined with the b = 
spectra of the foreground galaxies in each pair. These are used 
together to measure each foreground galaxy's CGM along two 
independent lines of sight. Composite spectra stacked accord- 
ing to galactocentric impact parameter allow us to measure the 
dependence of the line strength of several ionic species on im- 
pact parameter b from 0-125 kpc (physical). These lines are 
strongly saturated, so their strength is determined by a combi- 
nation of the covering fraction and the velocity spread in the 
absorbing gas. The fall-off of Wo with b implies f c (r) oc r" 7 , 
with 0.2 < 7 < 0.6, with Lya and CIV having smaller val- 
ues of 7 compared to (e.g.) Sill and SilV. For each ob- 
served species, the strength of absorption declines much more 
rapidly beginning at b ~ 70 - 90 kpc, with the exception of 
Lya which remains strong to b ~ 250 kpc. Combining the 
known space density of the galaxies in our sample (i.e., in- 
cluding only those with apparent magnitude 1Z < 25.5) with 
the average cross-section for absorption accounts for ~ 45% 
of all intergalactic C IV absorption with W (1548) > 0.15 A 
and ~ 70% of LLSs at z ~ 2 observed along QSO sightlines. 

6. We have proposed a simple model constrained by a com- 
bination of Wo(b) from the galaxy pairs and the kinematics of 
IS species observed toward each galaxy's own stars (b = 0). 
The IS absorption line shapes depend on a combination of 
v ut(r), the radial dependence of outflow velocity which in 
turn depends on the cloud acceleration a(r) and f c (r, v), the 
radial dependence of the covering fraction, which depends on 
cloud geometry and is inferred from Wo{b). We show that a 
self-consistent model using parametrized forms of a(r) and 
f c (r) can simultaneously reproduce the observed line shapes 
and strengths (apparent optical depth vs. velocity) as well 
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as the observed Wo(b) relation from the galaxy-galaxy pair 
measurements. In the model, higher velocity gas is located 
at larger galactocentric radii, and gas clouds are accelerated 
to ~ 800 km s _I before the absorption strength (i.e., covering 
fraction) is geometrically diluted. There is little evidence in 
the line profiles for stalling or infalling wind material, and 
the transverse sightlines show metal-enriched gas to at least 
~ 125 kpc 29 Taken together, these suggest that high velocity 
outflows from ~ L* galaxies at z ~ 2 - 3 are able to deposit 
enriched material to very large radii. They also suggest that 
the observed metals are causally related to the observed galax- 
ies, rather than remnants of winds generated by other (earlier) 
galaxies deposited in the nearby volume. 

7. The observed v max ~ 800 km s" 1 that appears to be a gen- 
eral property of the galaxies in our spectroscopic sample ex- 
ceeds v esc (r V i r ) estimated for galaxy halo masses of ~ 9 x 10 11 
M Q , the average halo mass expected given the observed space 
density and clustering of the z — 2 galaxies. Similarly, the 
galactocentric radius at which the strength of low-ionization 
absorption line species begins to decrease rapidly is similar 
to the estimated virial radius of ~ 80 kpc for the same halo 
mass. 

8. We consider the observations in the context of "cold 
accretion" or "cold flows", which are believed to be active 
over the same range of galactocentric radii and redshifts, for 
galaxies in the same range of total mass as those in our sam- 
ple (4 x 10 11 ZM tot < 2 x 10 12 M Q ). There is possible ev- 
idence for the presence of infalling gas in the far-UV spec- 
tra of galaxies with greater than the median baryonic mass 
Mbar = 4 x 10 10 M©, while the lower-mass sub-sample shows 
no evidence for the expected redshifted components of IS ab- 
sorption lines. The kinematics of Lya emission also seem in- 
consistent with the expectations for accreting material, while 
they are as expected in the context of the picture of out- 
flows we have presented. We expect that the same gas seen 
in absorption to b ~ 50- 100 kpc also acts as a scattering 
medium for escaping Lya photons initially produced in the 
galaxy H II regions. These Lya "halos" are probably a generic 
property of all high redshift LBG-like galaxies, the more ex- 
treme of which have probably already been obser ved as "Lya 
Blobs" (( Steid el et al .11200 0: Matsu da et alll2004l) '). More typ- 
ical LBGs likely exhibit this diffuse Lya emission, which is 
of low surface brightness and therefore difficult to observe 
without extremely deep Lya images. 

9. The inferred mass of cool gas in the CGM (within ~ 100 
kpc) is comparable to the sum of cold gas and stars in the 
inner few kpc; together they account for £ 50% of the total 
baryonic mass (~ 1.5 x 10 11 M Q ) associated with the aver- 
age M~9x 10 11 M Q dark matter halo. If the qualitative 
picture of the CGM we have proposed is correct, the wind 
material would still be traveling at v out > v e , c (r V! >) when it 
crosses the virial radius at r ~ 80 kpc. Even if v ou , ~ v esc , 
the low-ionization wind material would reach r ~ 100 kpc 
within ^150 Myr of the onset of the star formation episode 
if it is slowed only by gravity and does not accumulate a large 
amount of swept-up ISM on its way. Since we have shown 
that the CGM is outflowing and itself comprises a large bary- 
onic reservoir, it is possible that in most directions there is 
little to impede outflows from propagating ballistically into 



the IGM. 

The use of resolved background sources for studying the 
CGM in absorption is qualitatively different from using QSO 
sightlines. It seems clear that the multi-phase ISM observable 
in a number of ionic species indicates differential changes 
in covering fraction of the outflowing gas as a function of 
ionization level. This means that the radial behavior of the 
strength of saturated transitions is providing information on 
the "phase-space" density of gas having physical conditions 
amenable to the presence of each ion. Because absorption 
seen in the spectra of background galaxies averages over a 
~ 2- 3 kpc region, as compared to background QSOs, whose 
apparent angular sizes are ~ 5 orders of magnitude smaller, 
statistics such as IS line strength, covering fraction, and ve- 
locity range vs. impact parameter, are much less subject to 
large variance along a single sightline. 

All of the galaxy spectra used in this paper were obtained 
as part of large surveys, and therefore have low to moderate 
spectral resolution and S/N. It is possible with current gen- 
eration telescopes and instruments to obtain s pectra of rela- 
tively high S/N and moderate resolution (e.g., Hhapley et al. 
2006) in 10-20 hours integration time (i.e., ~ 10 times longer 
than most in the current samples). However, spectra similar 
to those currently possible only for strongly lensed galaxies 
(R ~ 5000, S/N > 30) will be routine at TZ ~ 24-24.5 us- 
ing future 30m-class telesco pes equipped with state of the art 
multi-object spectrometers (Steid el et al.l 120091) . When such 
capability is realized, full observational access to the three- 
dimensional distribution of gas near to and between galax- 
ies will revolutionize the study of baryonic processes in the 
context of galaxy formation. The results described above are 
intended as an illustration of the application of rest-far-UV 
galaxy spectra toward a better understanding of feedback pro- 
cesses and the CGM/IGM interface with the sites of galaxy 
formation. The feasibility of sharpening our understanding 
and interpretation with higher quality data in the future is very 
encouraging. 
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29 | Adelberger et d] j200l) and^ delberser et al. 1 2005b) have argued that at approximately this distance due t0 a number of separate considerations, 

there is evidence that the metals affect regions or ~ 300 kpc (proper) based on 
C IV-galaxy cross-correlation, and that outflows would be expected to stall 
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